ISSN 1175-5326 (print edition) 



Zootaxa 3741 (4): 401-458 
www.mapress.com/zootaxa/ 


Copyright © 2013 Magnolia Press 


Article ZOOTAXA 

ISSN 1175-5334 (online edition) 


http://dx.doi.Org/10.11646/zootaxa.3741.4.l 

http://zoobank.Org/urn:lsid:zoobank.org:pub:0D2276A9-5295-4622-BCE8-A7D996A30367 


Using various lines of evidence to identify Chironomus species 
(Diptera: Chironomidae) in eastern Canadian lakes 


ISABELLE PROULX 1 , JON MARTIN 2 , MELISSA CAREW 3 & LANDIS HARE 14 

1 Institut national de la recherche scientifique - Centre Eau Terre Environnement, 490 rue de la Couronne, Quebec City, Quebec, G1K 
9A9, Canada. E-mails: isabelle_proulx@hotmail.com & landis@ete.inrs.ca 

2 Department of Genetics, University of Melbourne, Melbourne, Victoria 3010, Australia. E-mail: j.martin@unimelb.edu.au 
3 Centre for Aquatic Pollution Identification and Management, Bio21 Institute, University of Melbourne, Melbourne, Victoria 3010, 
Australia. E-mail: mecarew@unirnelb.edu.ca 
4 Corresponding author 


Table of contents 


Abstract.401 

Introduction .402 

Methods. 403 

Results and discussion. 411 

Species delimitation and identification. 411 

Species descriptions and taxonomic status.421 

Chironomus (Chironomus) ‘tigris’ . 421 

Chironomus (Chironomus) staegeri Lundbeck (1898) .421 

Chironomus (Chironomus) frommeri Sublette and Sublette (1971).422 

Chironomus (Chironomus) cucini Webb (1969) .423 

Chironomus sp. NAII.423 

Chironomus sp. NAIII (possibly C. decumbens (Malloch 1934)). 424 

Chironomus sp. NAI (C. anthracinus- group). 424 

Chironomus (Chironomus) anthracinus Zetterstedt (1860).425 

Chironomus (Chironomus) ends Shobanov (1989) .425 

Chironomus (Chironomus) plumosus Linnaeus (1758).426 

Chironomus (Chironomus) maturus Johannsen (1908).427 

Chironomus (Chironomus) decorus- group sp. 2 Butler et al. (1995).427 

Chironomus (Chironomus) harpi Sublette (in Wulker et al. 1991) .428 

Chironomus (Chironomus) bifurcatus Wuelker, Martin, Kiknadze, Sublette and Michiels (2009).428 

Chironomus (Chironomus) dilutus Shobanov, Kiknadze and Butler (1999).429 

Chironomus (Chaetolabis) ochreatus Townes (1945) .430 

Chironomus (Chaetolabis) nr. atroviridis (sp. 2i) Martin (2013).430 

Morphological key to larvae of the Chironomus species collected in our study. 431 

Approaches used to delimit Chironomus species. 432 

Overall conclusions and recommendations for identifying Chironomus species.434 

Acknowledgements .434 

References .435 

Supplementary information (Tables S1-S5).440 


Abstract 

Chironomus Meigen (Diptera, Chironomidae) larvae are usually the largest sediment-burrowing chironomids, and as such 
often constitute a major part of the freshwater infaunal biomass. However, use of this genus in ecological, environmental 
and paleoecological studies is hampered by the fact that Chironomus larvae are difficult to identify to species because the 
larvae of many species are morphologically similar. We used a combination of morphological, cytological and genetic 
techniques to distinguish Chironomus larvae collected from 31 water bodies located in eastern Canada, producing 17 dis¬ 
tinguishable groupings. These groups of larvae were ultimately identified as belonging to 14 known species (C. anthraci- 
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mis, C. bifurcatus, C. cucini , C. decorus- group sp. 2, C. dilutus, C. ends, C. frommeri, C. harpi, C. maturus, C. nr. 
atroviridis (sp. 2i), C. ochreatus, C. plumosus, C. staegeri and C. ‘tigris ’) and three other species that remain unidentified 
(C. sp. NAI-III). No single approach served to delimit and identify larvae of all 17 Chironomus species that we collected. 
Although we expected that morphological criteria alone would be insufficient, our results suggest that DNA barcoding, 
using either the mitochondrial coxl or the nuclear gb2p gene, was also inadequate for separating some Chironomus spe¬ 
cies. Thus we suggest that multiple approaches will often be needed to correctly identify Chironomus larvae to species. 

Key words: Chironomus, morphology, cytology, DNA barcoding, coxl, gb2p, Canada 


Introduction 

The insect genus Chironomus Meigen (Diptera, Chironomidae) is found in fresh waters on all continents except 
Antarctica. It includes several hundred species, now classified into three subgenera ( Chaetolabis , Chironomus, 
Lobochironomus ) (the subgenus Camptochironomus is no longer recognized—see blether (2012)), as well as other 
species that are yet to be described (Ryser et al. 1985; Ashe & Cranston 1990; Martin 2013). In lakes from the 
tropics (Hare & Carter 1986), to the temperate (Jonasson 1972), to the Arctic (Butler 1982), Chironomus larvae are 
usually the largest sediment-burrowing chironomid and often represent a major part of the infaunal biomass. Thus 
Chironomus larvae can be an important source of food for fish and are widely used in ecological (Jonasson 1972), 
environmental (Martin et al. 2008) and paleoecological (Brooks et al. 2007) studies of fresh waters. If we are to 
understand their roles in aquatic ecosystems, it is important to be able to correctly identify Chironomus species. 

The identification of Chironomus larvae to species can be problematic because there are few conspicuous 
morphological differences among many Chironomus species (Lindeberg & Wiederholm 1979). As a result, larvae 
are often referred to simply as Chironomus spp. (Nyman et al. 2005) or at best are grouped into types according to 
the presence and form of their abdominal tubules (Shobanov et al. 1996) or the shape of their mouth parts (Brooks 
et al. 2007). Such groupings can limit the use of Chironomus larvae in ecological, environmental and 
paleoecological studies because behavioural and ecological differences among species are often important. For 
example, cadmium concentrations in sympatric Chironomus species can vary by an order of magnitude because of 
differences in their feeding habits and consequent contaminant exposure (Martin et al. 2008; Proulx & Hare 2008, 
2013). Pooling such species would clearly limit their use as contaminant biomonitors. If we cannot correctly 
identify Chironomus larvae to species, then it is difficult to use them to infer environmental impacts. 

In early studies, features of the head capsule and abdominal tubules were used to identify Chironomus larvae to 
species (Johannsen 1937). Subsequently, Chironomus species were also separated on the basis of the structure of 
polytene chromosomes located in their salivary glands (Keyl 1962; Martin 1979; Wulker et al. 1989). In the last 
decade or so, genetic techniques have been used to supplement these earlier taxonomic methods. For example, the 
polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) approach has been used to 
produce species-specific DNA profiles that can differentiate among Chironomus species (Carew et al. 2003; 
Sharley et al. 2004). This technique involves first amplifying specific genes or regions with PCR, and then 
digesting the resulting PCR amplicons with restriction endonucleases. Restriction endonucleases cut PCR 
amplicons differentially based on nucleotide differences in their DNA sequence, thereby generating a species- 
specific RFLP or DNA profile. The DNA profiles are visualised by gel electrophoresis as DNA fragments of 
different lengths. Although this method is inexpensive and useful for screening large numbers of individuals, it 
only examines a subset of the variation present in PCR amplicons (Pfrender et al. 2010). 

Another genetic technique used to separate and identify species is DNA sequencing (also known as Sanger 
sequencing) of PCR amplicons. This technique, referred to as DNA barcoding when used for identifying species, is 
more exact than PCR-RFLP as it detects all nucleotide differences. The standard gene used for DNA barcoding is 
the 3’ end of the mitochondrial cytochrome oxidase subunit I {coxl', Hebert et al. 2003). Advantages of using the 
coxl gene are that universal primers are able to amplify this gene from many animal groups (Folmer et al. 1994) 
and sequence variations in coxl can be used to discriminate among many closely-related species (Hebert et al. 
2004a). In insects, DNA barcoding using the coxl gene has been used to identify species from a range of groups 
including the Collembola (Hogg & Hebert 2004), the Ephemeroptera (Ball et al. 2005; Elderkin et al. 2012), the 
Coleoptera (Davis et al. 2011)andthe Chironomidae (Carew et al. 2007; Ekrem et al. 2007; Pfenninger et al. 2007; 
Sinclair & Gresens 2008; Ekrem et al. 2010; Carew et al. 2011; Stur & Ekrem 2011). Although coxl sequences can 
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be used to separate the majority of species, its mitochondrial origin is problematic for Chironomus because some 
species are known to hybridize (Martin 2011). Therefore, including sequence data from additional nuclear markers 
whose mode of inheritance differs from mitochondrial genes is required (Guryev et al. 2001; Martin et al. 2002; 
Martin 2011). To this end, the nuclear gene globin 2ft ( gb2(J) has been used in several studies on Chironomus 
species (Kao et al. 1994; Hankelneta/. 1997; Guryev et al. 2001; Guryev & Blinov 2002; Martin et al. 2002). 

We applied morphological, cytological and genetic techniques to identity Chironomus larvae collected in 31 
water bodies in eastern Canada to determine what combination of techniques would allow us to accurately identify 
the Chironomus species. To date, very few studies have used multiple techniques to discriminate among 
Chironomus species. We anticipate that the results of our study will be useful to those wishing to identify North 
American Chironomus species and will provide useful tools to those wishing to identify Chironomus species on 
other continents. The ability to accurately identify Chironomus species should facilitate future ecological and 
environmental studies in this and other geographical zones. 


Methods 

Collection and dissection of larval Chironomus. We collected fourth-instar Chironomus larvae from 31 water 
bodies (Table 1) located in the provinces of Quebec (near Quebec City, Rouyn-Noranda and Trois-Rivieres) and 
Ontario (near Sudbury), Canada. The collection period extended from ice-off in late spring (May) to early summer 
(June) in various years from 2006 to 2011 (Table 1). Exceptionally, in Lake Bedard, Chironomus were collected at 
the end of the summer (September). Sediments were collected using an Ekman grab and sieved through a net to 
eliminate fine sediment and retain Chironomus larvae. Larvae were preserved in 94% ethanol. 

Depth and water chemistry were measured at each collecting site (Table 1). Water samples were filtered in situ 
using diffusion samplers (Ponton & Elare 2009). Dissolved organic carbon (DOC) was measured by combustion 
and transformation into C0 2 (TOC-VCPEI, Shimadzu, Columbia, MD, USA) and magnesium (Mg) and calcium 
(Ca) concentrations were measured by inductively coupled plasma-atomic emission spectrometry (ICP-AES, Vista 
AX CCD, Varian, Mississauga, Ontario, Canada). Quality assurance of water chemistry measurements was assured 
through the use of blanks and appropriate standard reference materials. At the time of sampling, the water columns 
of all lakes were well mixed and oxygenated. 

Plead capsule and terminal abdominal segments with attached tubules were separated and kept for 
morphological studies, whereas the rest of the body was retained for genetic analyses. In addition, three individuals 
of each species (as determined by genetic analyses and morphology) were chosen at random for examination of 
their polytene chromosomes. For this purpose, the thoracic segments containing the salivary glands were preserved 
in a 3:1 mixture of 94% ethanol to glacial acetic acid (not all specimens showed chromosomal patterns of sufficient 
quality for species identification). Exceptionally, 41 larvae of C. entis and C. plumosus were examined 
cytologically to validate genetic results for these species. 

Genetic analyses. Polymerase chain reaction—restriction fragment length polymorphism (PCR-RFLP) 
analysis was performed on Chironomus larvae collected in 2006 and 2007. Specimens of each PCR-RFLP profile 
were subsequently sequenced for further DNA analysis. Specimens collected in subsequent years were sequenced 
directly. 

DNA extraction. DNA was extracted from larvae using the modified Chelex method (Carew et al. 2003). 
Briefly, the larval body (minus the head and terminal segments) was dried using a paper towel and placed in a 0.5 
mL plastic microcentrifuge tube. Individuals with large amounts of sediment in their gut were avoided as this can 
inhibit the PCR (Carew et al. 2003). Tubes were immersed in liquid nitrogen and the contents crushed into a 
powder using a pestle, 400pL of suspended 5% Chelex-100 resin (BioRad) was added and samples were incubated 
at 90 °C for 30 min. Extracts were stored at -20 °C until required for the PCR procedure. 

Polymerase chain reactions. PCR amplification of portions of the coxl and the gb2p genes was carried out in 
a 40 pL reaction mixture containing: lx PCR pH 8.8 buffer (20 mM Tris-HCl, 10 mM (NH 4 ) 2 S0 4 , 10 mM KC1, 2 
mM MgS0 4 , 0.1% Triton X-100; New England Biolabs (NEB)), 200 pM each of deoxynucleotide triphosphate 
(dNTPs), 0.4 mg/mL of bovine serum albumin (BSA), 0.5 pM of forward and reverse primers (see Table 2), 1 unit 
of Taq DNA polymerase (NEB), and 5 pL of Chelex DNA extraction supernatant taken from just above the resin 
after centrifugation at 15,000 relative centrifugal force for 2 min (Carew et al. 2003). All gh2fi gene primers tested 
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TABLE 1. Location, year and depth of collection, as well as water chemistry and trophic status of the water bodies studied. 
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in our study are listed in Table 2. The primers used to amplify the gb2(3 gene for each species are given in Table 4. 
For the coxl gene, the PCR thermal regime consisted of an initial denaturation cycle of 94 °C for 3 min; followed 
by 40 cycles of denaturation at 94 °C for 15 s, annealing at 45 °C for 45 s, elongation at 72 °C for 1 min; and one 
cycle at 72 °C for 1 min. The PCR thermal regime for the gb2/ J > gene was the same, but in some instances an 
annealing temperature of 50 °C was used. All PCRs had a negative control with no DNA template added. PCR 
products were verified by electrophoresis on a 1.5% Tris-Acetate-EDTA agarose gel with ethidium bromide. PCR 
product sizes were estimated using Flyper Ladder II (Bioline). 


TABLE 2. Primers used in this study. 


Gene 

Primer 

(forward (for) or reverse (rev)) 

Sequence (5’-3’) 

Reference 

coxl 

911 (for) 

912 (rev) 

TTT CTAC A A AT C ATA A AG ATATT GG 

TAAACTTCAGGGTGACCAAAAAATCA 

Folmer et al. (1994) 

gb2f 

wyklb (for) 

GAYATCCTTTACTACTYTT 

Modified version of Kao et al. 
(1994) wykl primer 


wyk4 (rev) 

wyk3 (rev) 

GACCTTGTGTCCAGGC 

GTGTTTCCATAGCTGGC 

Kao et al. (1994) 


2/?-B (for) 

2f -A (rev) 

2/?con for (for) 

2/?con rev (rev) 

GATATCCTTTACTACATC 

CG AT GT C A ATA A ATAC ATG 

CCAGACATCATGGCTAA 

CTTGACAACATCTTCGAC 

Hankelneta/. (1997) 


Polymerase chain reaction-restriction fragment length polymorphism analysis. Digest enzymes were 
chosen based on previous publications on chironomid identification using PCR-RFLP (Carew et al. 2003; Sharley 
et al. 2004; Carew et al. 2005; Carew et al. 2007). The PCR products from coxl were cleaved using restriction 
endonucleases with 4 base pairs (bp) (Alu I, Rsa I, Taq I) and 6 bp (Hha I, Hinf I, Ssp I) recognition sites. All 
digests were carried out as described by Carew et al. (2003) in a 20 pL reaction mixture containing: 10 pL of PCR 
product, lx recommended buffer, 0.1 mg/mL BSA and variable units of restriction endonucleases (3 units for Alu I, 
Rsa I and Ssp I, 4 units for Taq I and 6 units for Hha I and Hinf I; NEB). Restriction digests were incubated at 37 °C 
overnight, with the exception of Taq I, which was incubated at 65 °C for 3h. Digest products were separated via 
electrophoresis for 2h at 100V on a 3% agarose gel stained with ethidium bromide and observed under UV light. 
The size of digest fragments was estimated with a 50 bp ladder (Promega). Fragment sizes below 100 bp were 
ignored, as they were not always clearly discernible on the agarose gels. To verify results obtained from these 
digests, we simulated digests of the corresponding coxl sequences using the New England BioLabs NEBCutter 
V2.0 program (http://tools.neb.com/NEBcutter2/). Simulation digests were also performed on the coxl sequences 
of larvae for which RFLP digests were not made. 

DNA sequencing analysis. The coxl (709 bp) and gb2/3 (332-394 bp) gene products were purified and 
sequenced in both directions using the forward and reverse primers used for PCR amplification by Macrogen 
(Seoul, Korea) or by the research center at the Centre hospitalier universitaire de Quebec (Quebec, Canada) on an 
ABI3730 XL automatic DNA sequencer (Applied Biosystems) and were aligned using BioEdit 7.1.3.0 (Flail 1999). 
All sequences used for DNA analyses were submitted to GenBank (KF278208-KF278447; KF278449-KF278450). 
Sequences from coxl were aligned using CLUSTAL W (Thompson et al. 1994). Due to the presence of introns in 
some species, gb2f -sequences were aligned manually according to Flankeln et al. (1997). Sequences were 
analyzed in MEGA 5.05 (Tamura et al. 2011). Primer sequences for each gene were excluded from the analysis. 
Since our goal was to separate Chironomus species based on sequence similarities, rather than to infer interspecific 
phylogenetic relationships, identification trees (ID-trees) based on coxl -sequences and g/G/Fsequences were built 
using the Neighbor-Joining (NJ) (Saitou & Nei 1987) algorithm. The pairwise distances were calculated from the 
Kimura 2-parameter (K2P) model (Kimura 1980), which is best suited when distances are low (Nei & Kumar 
2000), as in our study. Bootstrap analysis was performed with 1000 replicates. For the coxl identification (ID) tree, 
Polypedilum aviceps. Drosophila affinis and Glyptotendipes lobiferus sequences from GenBank were added as 
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outgroups. In the case of th egb2ft ID-tree, because this gene is quite variable and because the only really conserved 
regions are also conserved in the globin genes 7A ( gb7A ) and 9 ( gb9 ) (Hankeln et al. 1997), Chironomus species 
sequences of the gb7A and gb9 genes were also added to make sure that all sequences obtained for our specimens 
were of the gb2[ J > gene. The gb2fl primers did in fact amplify the gb7A gene from C. ( Chaetolabis ) nr. atroviridis 
(sp. 2i). Pairwise intraspecific and interspecific nucleotide-sequence divergences were also calculated for all 
sequences using the K2P model in MEGA 5.05. Since some authors have expressed reservations about using 
divergence thresholds to separate species (DeSalle et al. 2005), including those of Chironomus (Martin 2011), we 
also used specific base differences to quantify differences between some closely-related species. 

Morphological analysis. Larval length was measured under a dissecting microscope and head capsule width 
(at the level of the eyes) and abdominal tubule lengths were measured using a microscope linked to an image- 
analysis system. 

Chironomus were sorted according to larval type on the basis of the presence or absence, types and length of 
abdominal tubules using a dissecting microscope. Although there have been several attempts to classify 
Chironomus larvae into types based on the morphology of their tubules (Harnisch 1942; Andersen 1949; Lindeberg 
& Wiederholm 1979; Shobanov et al. 1996; Shobanov 2002), the definitions in these schemes have been 
inconsistent and in some cases contradictory. In general we have returned to the original scheme of Harnisch 
(1942), for which many of the types were well-illustrated by Andersen (1949). To this scheme we have introduced 
two additional types, bathophilus-type and melanotus-type, from the more recent scheme of Shobanov (2002). 
While the early schemes clearly recognised the coiled nature of the ventral tubules of C. thummi (now C. riparius ) 
and C. plumosus, this aspect was lost in later classifications, such as those of Lindeberg and Wiederholm (1979) 
and Shobanov (2002), which were based only on the length of the tubules. We have found the distinction between 
coiled and relatively straight tubules to be consistent within species, and also a useful distinction among species. 
The bathophilus- and melanotus-types can fill this gap but only if the Shobanov (2002) definition is broadened to 
cover larvae with long ventral tubules, but without the typical coiling of those seen in thummi-type and plumosus- 
type. It should be noted that Shobanov (2002) introduced the melanotus-type to replace the anthracinus-type 
because C. anthracinus does not have an anthracinus-type larva (i.e. with lateral tubules), but is a typical thummi- 
type (i.e. without lateral tubules) (see larval description in the “Results and discussion” section). Further, the 
anthracinus-type, along with the semi-thummi- or semi-bathophilus-types, were intended to define in part larvae 
with very small lateral tubules, but we have found it difficult to draw a clear line between short and long lateral 
tubules and so have ignored this criterion. Our amended version of the larval classification is presented in Table 3. 


TABLE 3. Classification of Chironomus larval types. See Fig. 1 for illustrations of ventral tubules. 


Larval type 

Pair of lateral 

Two pairs of ventral tubules on IP 1 ’ segment 3 


tubules on 10"’ 
segment 

Anterior pair 

Posterior pair 

salinarius 

absent 

absent 

absent 

halophilus 

absent 

absent or short 

short 

bathophilus 

absent 

straight; long 

straight; long 

fluviatilis b 

absent 

slightly curved, coining to a point at ends; 

slightly curved, coming to a point at ends; 

thummi 

absent 

long 

with elbow; long 

long 

coiled; long 

reductus 

present 

absent 

absent 

semireductus 

present 

straight; short 

straight or may be slightly curved; short 

melanotus 

present 

straight or slightly curved; long 

straight or slightly curved; long 

plumosus 

present 

with elbow; long 

coiled; long 

a long: ventral tubules > the width of 11 th segment 

short: ventral tubules < the width of 11 th segment 
b Often hard to distinguish from bathophilus-type 
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FIGURE 1. Ventral tubules of the various larval types: salinarius and reductus (a), bathophilus and melanotus (b), fluviatilis 
(c), thummi and plumosus (d), as well as semireductus (e). 



I 


Pale 


Dark 


FIGURE 2. Dorsal view of larval head-capsules showing variation among species in the color of the frontoclypeus from pale 
(C. sp. NA1II) to dark (from left to right: C. ‘tigris C. cucini, C. dilutus). 
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p a j e Slightly Posteriorly Strongly Completely 

darkened darkened darkened darkened 

FIGURE 3. Ventral view of larval head-capsules showing variation among species in the color of the gula from pale (C. 
maturus) to slightly darkened (C. harpi), to posteriorly darkened (C. nr. atroviridis (sp. 2i)), to strongly darkened (C. staegeri ), 
to completely darkened (C. ‘tigris’). 



Type A 

Type B 

TypeC 

Type D 

Outer teeth almost 

Outer teeth only 

Outer teeth almost 

Outer teeth completely separated from 

completely fused to 

partially separated 

completely 

middle tooth. 

middle tooth 

from middle tooth. 

separated from 
middle tooth. 



FIGURE 4. Types of mentum middle trifid tooth: (type A) C. nr. atroviridis (sp. 2i), (type B) C. cucini, (type C) C. staegeri, 
(type D) C. maturus (left) and C. plumosus (right). 



Type I 

4 th lateral tooth the same size as 
the 3 rd and 5 th lateral teeth 


Type II 

4 th lateral tooth reduced, it is 
the same height as the 5 th 
lateral tooth 



Type III 

4 th lateral tooth greatly reduced, 
it is shorter than the 5 th lateral 
tooth 


FIGURE 5. Types of mentum 4“’ lateral teeth (Webb & Scholl 1985; Vallenduuk & Moller Pillot 1997): (type I) C. plumosus, 
(type II) C. staegeri, (type Ill) C. cucini. 


Head capsules were separated into parts and mounted in Canada Balsam so as to determine the coloration of 
the frontoclypeus (Fig. 2) and the gula (Fig. 3) as well as the structure of: the central trifid tooth (Fig. 4) and 4 th 
lateral teeth of the mentum (Webb & Scholl 1985; Vallenduuk & Moller Pillot 1997) (Fig. 5), the mandibles (Fig. 
6), the pecten epipharyngis (Fig. 7) and the ventromental plates (Webb et al. 1985). We developed a classification 
scheme based on differences in the central trifid tooth of the mentum (Fig. 4) and teeth of the mandibles (Fig. 6) 
that is based in part on the previous classifications of Webb and Scholl (1985) and Vallenduuk and Moller Pillot 
(1997), but that better encompasses the range in variation we observed in these structures. For example, Webb and 
Scholl (1985) classified the central trifid tooth of the larval mentum according to the degree of fusion of its three 
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component teeth, the width of the middle tooth, and the height of the outer teeth relative to the middle tooth. We 
found that the latter two criteria varied substantially within species and so considered only the first of these three 
criteria for that character. We note that although the degree of sharpness of the teeth of the mentum has been used 
for separating some Chironomus species (Martin 2013), this feature varied widely within the species under study 
and thus we did not use it for separating our study species. Lastly, we used the coloration of the 3 rd inner 
mandibular tooth and its degree of fusion with the lower mandibular margin to classify larvae. 


Coloration of 3 rd inner tooth 


Pale 


> Partially dark —> Completely dark 


3 rd inner tooth 
partially to 
completely 
fused to the 
lower margin 



Type A 


3 rd inner tooth 
completely 
separated 
from the 
lower margin 



TypeD 




Type B 


Type C 


3 rd inner tooth 
partially dark 
and separated 
from lower 
margin 

(not found in our 
study) 

Type E Type F 



FIGURE 6. Mandible types as defined by the degree of darkening and separation of the 3 rd inner tooth: (type A) C. bifurcatus, 
(type B) C. sp. NAI (above) and C. ochreatus (below), (type C) C. nr. atroviridis (sp. 21), (type D) C. staegeri, (type F) C. 
plumosus. 


Cytological analysis. Isolated salivary glands were prepared for polytene chromosome analysis using the 
aceto-orcein method (Martin el al. 2006). Veronika Golygina (Institute of Cytology and Genetics, Novosibirsk, 
Russia) assisted in distinguishing cytogenetically between C. entis and C. plumosus. Preparations of these two 
species have been deposited at her Institute. Polytene chromosome mounts of the remaining species have been 
deposited, together with their respective head capsule mounts, at the Canadian National Collection of Insects, 
Arachnids and Nematodes in Ottawa, Ontario, Canada. 

Species delimitation and identification. Larvae were sorted according to their morphology, their coxl PCR- 
RFLP profiles (larvae collected in 2006 and 2007) and their coxl and gh2(l gene partial nucleotide sequences and 
then linked, via DNA sequences of cytologically-known species (either already in GenBank, or from karyotyping 
larvae also sequenced in this study) to recognized species. 
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Type A 

Pointed and rather 
uniform teeth 


Type B 

Pointed teeth with 
occasionally 
smaller teeth 
interspersed 
between larger 
teeth 


TypeC 

Rounded and rather 
uniform teeth 


Type D 

Rounded teeth with 
smaller teeth 
interspersed 
between larger 
teeth 


FIGURE 7. Types of teeth on the pecten epipharyngis: (type A) C. staegeri, (type B) C. anthracinus, (type C) C. dilutus, (type 
D) C. ochreatus. 


Results and discussion 

Species delimitation and identification 

Species identifications are performed on 4 th (final) instar larvae. We confirmed that Chironomus larvae were in the 
fourth instar by comparing the width of their head capsule (Table 6) to those of prepupal larvae and larval exuviae 
attached to pupae (data not shown). We did not measure the head capsule widths of Chironomus sp. NAII larvae, 
but these were undoubtedly 4 th instars because we collected them just prior to adult emergence. Fourth instar larvae 
can also be recognized by the presence of developing imaginal discs in the thorax and/or posterior abdominal 
segments (Wiilker & Gotz 1968; Ineichen et al. 1983). 

Analysis using PCR-RFLP of the coxl gene was performed on 296 larvae. The enzymes Ssp I, Hinf \, Rsa I and 
Taq I were used to cleave the partial coxl gene into different RFLP profiles (Table 4). These profiles were 
congruent with our groupings based on larval morphology (larval types and head-capsule features), with the 
exception of a single profile (Ssp I: 500,240; Hinf I: 710; Rsa I: 500,240; Taq I: 260,200,190) obtained for two 
larval types that differed in the coloration of their frontoclypeus. For these larvae, the partial coxl gene was cleaved 
with two additional restriction endonucleases, Hha I and Alu I, thereby creating three extra RFLP profiles. Results 
of these analyses are summarized in Table 4. To verify the accuracy of these results, larvae that included all of these 
RFLP profiles were sequenced. 

We sequenced the partial coxl gene of 59 larvae that included all 15 RFLP profiles, as well as that of 79 other 
larvae (Fig. 8). We also amplified and sequenced the partial gb2fi -gene of 83 larvae (Fig. 9). However, we were 
unsuccessful in obtaining the gb2[3 sequence for all Chironomus species (Table 4) despite modifying PCR 
conditions and testing several primer combinations (Table 2). 

For visualisation purposes only, one representative of each unique sequence was used to illustrate the 
relationship between species in the coxl (Fig. 8) and gb2/l ID-trees (Fig. 9). However, trees were also built using all 
sequences (including individuals that had identical gene sequences) which showed that using only unique 
sequences did not affect tree topology. Sequences were grouped into potential species according to molecular 
evidence (sequence clusters with bootstrap values >90% and sequence divergences of <4%) and larval morphology 
(see curly brackets in Figs. 8-9). Following this, species were identified through polytene chromosome analysis 
(71 larvae) and DNA barcoding. For cytological analyses, results of these identifications are given after the vertical 
line located to the right of the corresponding sequences in Figures 8 and 9. For DNA barcoding, Nearctic 
Chironomus coxl -sequences and g/>2/?-sequences from GenBank (http://www.ncbi.nlm.nih.gov) that grouped 
(bootstrap values >90% and sequence divergences <4%) with our sequences were added to our ID-trees. Moreover, 
we sequenced the coxl gene and/or the gb2/ J > gene of voucher specimens (Table S2) and added these sequences to 
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TABLE 4. Summary of molecular results for the coxl and gb2p genes for each Chironomus species. Coxl gene: cox7-RFLP fragment sizes in base pairs (bp); 710-bp 
fragments are uncut; number of individuals analysed from real digests and number of individuals who subsequently had their coxl gene sequenced. Gb2 gene: primers 
used to amplify and sequence the gb2f} gene; absence or presence of an intron (type I or type II) in the amplified gb2/3 sequence. ND; Not determined. 
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TABLE 5. Average (%) pairwise sequence divergence between species and within species (given diagonally in bold font) for the coxl and gb2 /? genes. For each 
species and each gene, the number of sequences analysed is given in parentheses. Interspecific sequence divergences within the intraspecific range of 
Chironomus species (coxl: 0-3% and gb2fi: 0-2%; see Table S3) are highlighted in grey. 
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C.staegeri 14.73 14.53 13.14 12.05 15.97 13.30 

C.'tigris’ 14.20 28.50 14.23 34.71 12.06 23.11 11.69 15.90 27.36 12.58 27.96 

’Sequence interspecific divergences range from 1 to 5%. 
















TABLE 5. (Continued) 
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TABLE 6. Larval morphology of fourth-instar larvae of the various Chironomus species collected in this study. 


Anterior 

margin of 

ventromental 

plates 

smooth 

smooth 

smooth 

smooth 

smooth 

relatively 

smooth 

smooth 

Pecten 

epipharyngis 

Mean 

no. of Type 

teeth (Fig. 7) 

(range) 

CQ < < 

1 > 1 ^ ^ 
^ (N ON CN 7 

(7 O’ - 1—1 Os 

< 

r© o (S' 

15 

(13- A 

20) 

CQ Q 

'O rn 00 '-O 4 (S' 

Mandible 

type 

3rd inner 
tooth 
(Fig. 6) 

A 

B 

A 

< 

< 

AorB 

C 

Mentum type 

Central 4th 

trifid lateral 

tooth teeth 

(Fig. 4) (Fig. 5) 

B III 

C II 

B II or III 

CQ 

B or C II 

B or C II 

A I 

Gula color 
(Fig. 3) 

posteriorly 

darkened 

posteriorly 

darkened 

posteriorly 
to strongly 
darkened 

posteriorly 

to strongly 

darkened 

strongly to 

completely 

darkened 

strongly to 

completely 

darkened 

darkened 

posteriorly 

Frontoclypeus 

color 

(fig- 2) 

pale or 
slightly 
darkened with 
lobed dark 
spot anteriorly 

dark 

longitudinal 
stripe with 
lobed dark 
spot medially 

pale 

pale 

pale 

pale 

pale 

Head width 
in mm 
(range) 

0.54 

(0.49-0.62) 

not measured 

0.48 

(0.43-0.54) 

0.49 

(0.42-0.59) 

0.56 

(0.51-0.60) 

0.64 

(0.60-0.67) 

0.57 

(0.56-0.57) 

Mean larval 
length in 
mm (range) 

18 

(15-21) 

12 

(11-13) 

14 

(10-18) 

13 

(11-16) 

14 

(11-16) 

18 

(13-21) 

19 

(17-20) 

Mean tubule length in mm (range) 

T . Anterior Posterior 

Lateral , , , , 

, . . ventral ventral 

u U es tubules tubules 

absent absent absent 

absent absent absent 

absent absent absent 

absent ^ ^ 

(0.9-1.9) (1.0-1.7) 

absent ^ 0.9 

0.18 (0.7-1.5) (0.5-1.4) 

(0.16-0.21) 

T) O 

P ~ OS ~ 

— p O 

0_ O 

(S' 0 

2 T 2 7 

— 1 p ^ p 

"S c 

0 a> 

c n co 

JO X> 

c3 c3 

Larval type 
(Table 3 and 
Fig. 1) 

salinarius 

salinarius 

salinarius 

bathophilus 

bathophilus, 

fluviatilis 

melanotus 

thummi 

thummi 

C 

26-31 

4 

47-52 

o 

1 

5-17 

22-27 

4 

Species 

C. cucini 

C. sp. NAII 

C. sp. NAIII 

C. bifurcatus 

C. decorus- 
group sp.2 

c. 

anthracinus 

C. nr. 
atroviridis 
(sp. 2i) 
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TABLE 6. (Continued) 
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our ID-trees (Figs. 8—9). C. acidophilus (Keyl 1960), C. calligraphus (Goeldi 1905), C. quinnitukqut (Martin et al. 
2010), C. sp. g (Martin 2013), C. sp. h (Martin 2013) and C. sp. u (Martin 2013) did not group with any of our 
collected species, so we could rule out these species as being any of our unrecognized species. Analysis of all the 
sequences together with simulation digests allowed us to identity the Chironomus species that had been previously 
separated using RFLP analysis. Simulation digests were performed on all of the coxl sequences obtained (see 
Table 4). Extra RFLP profiles were obtained from these simulation digests. 

On the basis of the genetic, morphological and cytological information that we obtained, we conclude that the 
404 Chironomus larvae that we collected represent 17 species, 14 of which were known while the status of three 
others remains uncertain. A detailed list of all the larvae analysed is presented in Table S1. In the following section, 
pertinent genetic, morphological and cytological information is presented for each of these species in the order that 
they are presented in Figure 8 (from top to bottom). Detailed genetic information, including cox /-RFLP sizes, the 
primers used to amplify the gb2p gene and whether or not the gh2f> gene was amplified is given in Table 4. The 
gb2[l gene of some Chironomus species includes an intron (type I or type II), whereas in others it is absent 
(Hankeln et al. 1997; Makarevich et al. 2000). This information is also given in Table 4. The intraspeciflc and 
interspecific divergences in coxl and gb2fi sequences are summarized in Table 5, with more detail presented in 
Table S3. Overall, intraspecific divergences for Chironomus species characterized by coxl range between 0 and 
3%, whereas intraspecific sequence divergences based on gb2fi ranged between 0 and 2%. A detailed 
morphological description for each species is found in Table 6. Where pertinent, we also discuss the status of each 
species collected and its relationship to other closely-related species that we did not collect in our study. Following 
our species descriptions, we present a morphological key to discriminate among larvae of the 17 Chironomus 
species that we collected. Species from our study area (Ontario and Quebec) that are not included in the key are 
given in Martin (2013). Lastly, the distribution and ecology of each of the Chironomus species that we collected are 
summarized in Table 7. 


TABLE 7. Lakes in which the various Chironomus species were collected as well as lake characteristics including 
sampling depth, pH and trophic status. Lake codes are given in Table 1. ND; not determined. 


Species 

Depth (m) 

pH 

Trophic status 

Water body 

C. anthracinus 

2-12 

4.5-8.5 

oligo to mesotrophic 

AR, HA, OS, PI, RA, RM, SI (2011) 

C. bifurcatus 

1.5-24 

2.7-7.8 

oligo to eutrophic 

AD, AR, DA, DP, KI, MC, OP, SJ, TI (2011) 

C. cucini 

9-35 

6.2-7.5 

oligo to mesotrophic 

BO, CL, OP, SJ, VA 

C. decorus- group sp.2 

1-5 

7.2-7.6 

oligo to mesotrophic 

AD, DF, DP, FO, OP, SC, SI (2011) 

C. dilutus 

1.5-5 

7.5-84 

eutrophic 

KE 

C. entis 

1-9 

7.1-8.3 

meso to eutrophic 

DA, DS, MN, OP, PE, 

C. frommeri 

1 

ND 

ND 

PO 

C. harpi 

1-4 

2.7-3.8 

oligotrophic 

AR 

C. maturus 

4-6 

ND 

ND 

BE 

C. nr. a troviridis (sp. 2i) 

1 

7.4 

meso to eutrophic 

MN 

C. ochre at us 

3 

7.7 

mesotrophic 

OP 

C. plumosus 

1-8 

6.8-8.5 

oligo to eutrophic 

AU, DA, DP, FO, KE, KI, MN, OS, PE, RO 

C. sp. NAI 

7.5 

6.8 

oligotrophic 

KA 

C. sp. NAI1 

4 

5.9 

oligotrophic 

SI (2007) 

C. sp.NAIII 

5-12 

7.1-7.9 

oligo to mesotrophic 

DA, HA, MC, RA, RM 

C. staegeri 

1-10 

5.9-8.0 

oligo to eutrophic 

CR, DA, DP, KA, KI, MC, OP, PO, SI 
(2007), SJ, TI (2007 and 2011) 

C. ‘tigris’ 

2-10 

5.9-8.0 

oligo to mesotrophic 

KA, MC, OP, SI (2007), SJ, TI (2007 and 
2011) 
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. tigris 


64r° plumosus (n=l / MC / KF278253) 
f- O plumosus (n=l /T1 / KF278252) 

O plumosus (n=IO / KA, MC, SI, SJ, TI / KF278242-KF27825I) | C ‘tigris ’ 

O plumosus (n= I / TI / KF278241) 

-A C. tigris '(K.F278350) 

—?2-| O plumosus (n=3 1 OP / KF278238-KF278240) | C. -tigris ' 

O plumosus (n=l / SJ / KF278268) 

- O plumosus (n-1 / MC / KF278254) 

A C.staegeri (KF278356) 

- O plumosus (n= 13 / CR, DA, DP, KI, MC, OP, PO, SI, SJ, TI / KF278255-KF278267) | C. stteger/ 
O plumosus (n=l / PO / KF278237) 

- O plumosus (n=2 / PO / KF278235-KF278236) | C. frommeri 
— Ac.sp. h(KF278349) 

O salinarius (n=l / CL / KF278290) | C. ci/c/W 
O salinarius (n=2 / CL / KF278288-KF278289) 

O salinarius (n=2 / CL, OP / KF278286-KF278287) 

O salinarius (n=2 / VA/ KF278284-KF278285) | C. cuctw 
pO salinarius (n=l / SI / KF278269) 

1-C 


C. staegeri 


C. frommeri 


C. cucirti 


O salinarius (n-1 / SI / KF278270) 

O salinarius (n=3 / FA, RA, RM / KF278281-KF278283) 
O salinarius (n=l / RM / KF278280) 

-O salinarius (n=4 / DA, RA / K.F278276-KF278279) 

- O salinarius (n=2 / HA, RA / KF278274-KF278275) 

O salinarius (n=3 / HA, RA / KF278271-KF278273) 

98J ° thumm i (n= 1 / KA / KF278219) 


C. sp. NAU 


C. sp. NAIII 


} 


NAI 


75 f" C 

*9|0 p| 

^Opi 


1 o Ihummi (n-1 / KA / KF278220) 

C. anihracinus-group (KF278344) 

O thummi (n=4/ HA, OS, PI / KF278221-KF278224) | C. anthracinus 
C. anthracinus (KF278343) 

O thummi (n=10 / AR, OS, PI, RA, RM, SI / KF278225-KF278234) | C. anthracinus 
— Ac. entis (KF278355) 
k C. plumosus (KF278354) 

99^ 5 JO semireductus/plumosus (n=2 / MN, PE / KF278216-KF278217) | C. plumosus 
emireductus/plumosus (n=l / MN / KF278218) | C. plumosus 
I C. enlis (AFI92195) 

— O semireductus/plumosus (n=l / AU / KF278214) | C plumosus 
— O semireductus/plumosus (n-1 / RO / KF278215) | C. plumosus 
„ IO semireductus/plumosus (n-2 / DA, OP / KF278213) | C. e/tf« and C. plumosus 
(emireductus/plumosus (n=l / PE / KF278212) | C. entis 
- O semireductus/plumosus (n=5 / DA, DP, DS, KI, OS / KF278208-KF278211) | C. entis ai 


} 


C. enlis and C. plumosus 


} 


gr.2 


C. maturus (DQ648204) 
plumosus (n=l / BE / KF278291) 
plumosus (n=I / BE / KF278292) 

B C. f/eeonu-group sp.2 (DQ648205) 

O fluviatilis (n-3 / DP, FO, SI / KF278299-KF27830I) 

O fluviatilis (n=l / AD / KF278298) 

O fluviatilis (n=3 / AD, OP / KF278295-KF278297) | C. rfecorrts-group sp.2 
O fluviatilis (n=l / SC / KF278294) 

L 0 fluviatilis (n=I / DF / KF278293) 

C. quinnitukqut (KF278347) 

991 — A C. harpi (KF278346) 

^ O plumosus (n=5 / AR / KF278302-KF278306) 

O bathophilus (n=l / DA / KF278315) 

O bathophilus <n=l 1 / AD, AR, DP, OP/ KF278316-KF278326) 

A C. bi/urcatus (KF278361) 

O bathophilus (n=l / SJ / KF278352) 

4T Ac. bifiuvatus (KF278345) 

•-Ac. bifurcatus (KF278353) 
l4jO bathophilus (n= I /KI / KF278307) 

8i y- O bathophilus (n=3 / SJ / KF278312-KF278314) 

TtjfO bathophilus (n=2/ MC / KF278310-KF278311) | C. bifurcatus 
O bathophilus (n=2 / MC, TI / KF278308-KF278309) | C. bifurcatus 
Ac. calligraphus (KF278357) 

Ac. sp. g (KF278348) 

A C. acidophilus (KF278358) 

C. dilutus (AF110162) 

O plumosus (n=2 / KE / KF278333-KF278334) 

C. pallidivittatus (AF110165) 

- A C. dilulus (KF278359) 

C. dilulus (AF110160) 

7 rO plumosus (n=2 / KE / KF278335-KF278336) 

5^r O plumosus (n=l / KE / KF278337) 

10 plumosus (n=4 / KE / KF278338-KF278341) 

76j A C. ( Chaetolabis) ochrealus (KF278351) 

O thummi (n=l / OP / KF278327) | C. ( Chaetolabis ) ochrealus 
O thummi (n-1 / OP / KF278328) | C. ( Chaetolabis) ochrealus 
A C. (Chaetolabis) nr. atroviridis (sp. 2i) (KF278342) 

C. (Chaetolabis) nr. atroviridis (sp. 2i) (KF278360) 

99 C O thummi (n= I / MN / KF278329) | C. (Chaetolabis) nr. atroviridis (sp. 2i) 

O thummi (n=3 / MN / KF278330-KF278332) | C. (Chaetolabis) nr. atroviridis (sp. 2i) - 

-B Pofypediium aviceps (JF287749) 

-B Drosophila qffinis (HQ979106) 

B Glyptotendipes lobiferus (DQ648213) 


C. cfecorws-group sp. 2 


C. harpi 


- C. bifurcatus 


C. dilutus and C. pallidivillalus 


}«■ 


C. nr. atroviridis (sp. 2i) 


FIGURE 8. Neighbor-joining identification tree (NJ ID-tree) based on partial coxl sequences and the K2P substitution model. 
Numbers on branches are bootstrap values >50%. o Sequences of Chironomus species collected from lakes in our study. Larval 
morpho-types are specified followed in parenthesis by: sample size (n = the number of individuals sequenced for each 
consensus sequence), lake abbreviations, and GenBank accession numbers. Some larvae were identified by examining their 
polytene chromosomes, and these results are indicated alongside the corresponding sequence next to the vertical line. ■ 
Sequences obtained from GenBank (species name and GenBank accession number in parenthesis). A Sequences obtained from 
cytologically identified reference Chironomus specimens (species name and GenBank accession number in parenthesis). 
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93 J A C. plumosus (KF278444) 

O semireductus/plumosus (n=2 / PE / KF278381-KF278382) | C. plumosus 
O semireductus/plumosus (n=1 /AU / KF278383) | C. plumosus 
| O semireductus/plumosus (n=5 / DA, MN, PE / KF278384-KF278388) | C. plumosus 
'O semireductus/plumosus (n=1 / AU / KF278389) | C. plumosus 
O semireductus/plumosus (n=1 / FO / KF278380) | C. plumosus 
O semireductus/plumosus (n=4 / AU. DP. FO. KE / KF278376-KF278379) | C. plumosus 
O semireductus/plumosus (n=1 / OS / KF278375) | C. plumosus 
O semireductus/plumosus (n=1 / PE / KF278374) | C. plumosus 
5 < -O semireductus/plumousus (n=10 / DA, DP, KE. OS, PE. RO / KF278362-KF278371) | C. plumosus 
O semireductus/plumosus (n=2 / KE / KF278372-KF278373) | C. plumosus 
O thummi (n=1 / AR / KF278426) 

O thummi (n=3 / KA / KF278390-KF278392) —> C. sp. NAI (type 11 intron) 

O thummi (n=1 / RA / KF278425) 

O thummi (n=1 / AR / KF278427) 

O thummi (n=1 / PI / KF278395) | C. anthracinus 
O thummi (n=1 / RA / KF278394) 

O thummi (n=1 / OS / KF278393) 

A c. ent/s (KF278445) 

semireductus/plumosus (n=7 / DA. DS, OP. PE / KF278410-KF278416) | C. entis 
C. entis (KF278446) 

O plumosus (n =6 / AR / KF278428-KF278433) —^ C. harpi (no intron) 


C. plumosus (type II intron) 


C. anthracinus (type II intron) 


" C. entis (type II intron) 


O bathophilus (n=2 / AR / KF278434-KF278435) 

O bathophilus (n=1 / AR / KF278436) 

O bathophilus (n=1 / DP / KF278437) 

,0 bathophilus (n=1 / OP / KF278438) 

J.O bathophilus (n=3 / DP, OP, Tl / KF278439-KF278441) | C. bifurcates 
lo bathophilus (n=2 / SJ / KF278442-KF278443) 

100 |O salinarius (n=2 1 VA/ KF278419-KF278420) | C. cucini 
Lf 


- C. bifiircatus (no intron) 


l O salinarius (n=2 / SJ I KF278417-KF278418) 

- A c. sp. u (KF278447) 

-O plumosus (n=4 / MC, OP Tl / KF278396-KF278399) 

_ 100 TO salinarius (n=2 / HA, RA/ KF278423-KF278424) 

' o salinarius (n=2 / RA, RM / KF278421-KF278422) 

-O plumosus (n=1 / BD / KF278400) 

-Ac. calligraphus (KF278449) 

,001 ■ C dilulus (AF110173) 

I C. dilutus (AF110174) 

IO plumosus (n=9 / KE / KF278401-KF278409) 

I C. pallidivittatus (AF110170) 

I C. pa»lV/«a/us(AF110167) 

I C. pallidivittatus (AF110169) 

B C. pallidivittatus (AF110168) 


} 


C. cucini (no intron) 


C. ‘tigris’ (type II intron) 

C. sp. NA1II (type II intron) 
C. maturus (no intron) 


} 


C. dilutus (no intron) 


- C. pallidivittatus (no intron) 


8j~C 


-O g67A - C. thummi (AJ003812) 

-A gb7A - C. nr. atroviridis (sp. 2i) (KF278450) 

-□ gb7A - C. metanotus (AJ003797) 

- □ gb9 - C. cingutatus (AJ003788) 

-□ gb9 - C. aprilinus (AJ003787) 

-□ gb9 - C. luridus (AJ003794) 

-□ gb9 - C. thummi piger (AJ003803) 


1_ 

99 

□ gb9 

r c 




□ gb9 - C. tentans (AJ003809) 


-□ gb9 - C. duplex (AJ003790) 

I □ gb9 - C. oppositus (AJ003800) 
1— □ gb9 - C. jacksoni (AJ003793) 


o.os 

FIGURE 9. Neighbor-joining identification tree (NJ ID-tree) based on partial gb2P sequences and the K2P substitution model, 
o Sequences of Chironomus species collected from lakes in our study. Larval morpho-types are specified followed in 
parentheses by: sample size (n = the number of individuals sequenced for each consensus sequence), lake abbreviations and 
GenBank accession numbers. Some larvae were identified by examining their polytene chromosomes, and these results are 
indicated alongside the corresponding sequence next to the vertical line. ■ Sequences obtained from GenBank (species name 
and GenBank accession number in parenthesis). ▲ Sequences obtained from cytologically identified reference Chironomus 
specimens (species name and GenBank accession number in parenthesis). □ Chironomus gb7A and gb) sequences obtained 
from GenBank were also added as outgroups (globin name, species name and GenBank accession number in parenthesis). 
CsChironomus gb7A sequence obtained from cytologically identified reference specimens (globin name, species name and 
GenBank accession number in parenthesis). 
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Species descriptions and taxonomic status 
Chironomus ( Chironomus ) Tigris’ 

(nomen nudum in Martin et al. (2008), for species C. sp. Ami of Kiknadze et al. (1993)). 

Material examined (Table SI): 66 larvae from Kasten Lake, McFarlane Lake, Silver Lake and Tilton Lake in 
Ontario as well as from Lake Opasatica and Lake St. Joseph in Quebec. 

DNA (Figs. 8-9 and Tables 4-5, S3). We were able to obtain PCR products for both the coxl and gb2f genes. 
Compared to the other Chironomus species (except C. sp. NAIII), the gb2f gene of C. tigris is 3 codons short 
immediately after the end of the 2 nd intron. C. ‘ tigris ’ sequences form distinct clades in both the coxl and gb2b ID- 
trees. Consequently either gene can be used to accurately separate and identify C. ‘tigris ’. However, the range of 
interspecific divergence between the coxl -sequences of C. ‘ tigris ” and C. staegeri (1-2%), as well as between C. 
‘tigris ’ and C. frommeri (2-3%), are within the intraspecific sequence divergence range of collected and reference 
Chironomus species (0-3%). Therefore, coxl sequence divergence values alone cannot be used to reliably separate 
C. ‘tigris’ from C. staegeri or C. frommeri. For the gb2b gene, we could not assess the interspecific sequence 
divergence between C. ‘tigris ’, C. staegeri and C. frommeri because we were unsuccessful in amplifying the gh2fl 
gene for C. staegeri and C. frommeri. In the coxl ID-tree, sequences of collected larvae cluster with the reference 
sequence of C. 'tigris ’, thus confirming the identification of this species. Chosen restriction enzymes for the coxl 
PCR-RFLP analysis correctly separated C. ‘tigris ’ larvae from the other Chironomus species. 

Morphology (Table 6). Our specimens are large sized plumosus-type larvae having their anterior ventral 
tubules longer than the posterior ventral tubules. The frontoclypeus of C. ‘tigris’ is dark-colored, which 
distinguishes it from larvae of the otherwise morphologically-similar C. staegeri and C. frommeri, which both have 
a pale frontoclypeus. Exceptionally, in larvae from some other regions, the frontoclypeus of C. staegeri is reported 
to be slightly darkened (Martin 2013). However, this criterion could prove to be less clear cut and other 
morphological features such as the length of the lateral tubules and the anterior margin of ventromental plates 
could be used to separate these species. 

Cytology. The cytology of the two larvae analyzed clearly indicates that this species is C. ‘tigris ’ since it is one 
of only two Chironomus species known to possess two polytene chromosomes. In this respect, it is clearly distinct 
from C. staegeri and C. frommeri which possess three and four chromosomes, respectively. The arm combination 
of C. ‘tigris ’ chromosomes is GAB, FEDC and its chromosomes are described in Martin et al. (1974), Butler et al. 
(1995, C. sp. r), Kiknadze et al. (1993, C. sp. Ami) and Martin (2013). 

Distribution and ecology (Table 7). This species has been previously reported from lakes in Minnesota, 
Ontario, Quebec and Wisconsin (Butler et al. 1995; Martin et al. 2008; Martin 2013). We found C. ‘tigris’ in 
oligotrophic to mesotrophic lakes of pH 5.9-8.0. At all sites where C. ‘tigris ’ was collected, C. staegeri was also 
present. However, the reverse was not necessarily the case. Lakes in which C. staegeri was present and C. ‘tigris ’ 
was absent tended to be eutrophic (with the exception of Crooked Lake), which suggests that C. ‘tigris ’ larvae 
prefer less productive systems. Their northerly distribution in North America may reflect this fact. We collected C. 
‘tigris ’ at water depths varying from 2-10 m, although it can live at greater depths (20 m; Butler et al. 1995). 


Chironomus ( Chironomus) staegeri Lundbeck (1898) 

Material examined (Table SI): 44 larvae from Lake D’Alembert, Lake Duprat, Lake Kinojevis, Lake Opasatica, 
Lake St. Joseph and an unnamed pond in Quebec as well as from Crooked Lake, McFarlane Lake, Kasten Lake, 
Silver Lake, and Tilton Lake in Ontario. 

DNA (Figs. 8-9 and Tables 4-5, S3). In the coxl ID-tree, C. staegeri sequences form a distinct clade. 
Sequences of collected larvae cluster with the reference sequence of C. staegeri, thereby confirming the 
identification of this species. We were not successful in amplifying the gb2f gene for this species. Coxl 
interspecific sequence divergences between C. staegeri and C. ‘tigris’ (1-2%) as well as C. staegeri and C. 
frommeri (3-4%) are within the intraspecific divergence range of Chironomus species assessed in this study (0- 
3%). Therefore, coxl sequence divergence values cannot be used to separate C. staegeri from C. ‘tigris' or from C. 
frommeri. Restriction enzymes for the coxl PCR-RFLP analysis correctly separated C. staegeri from the other 
Chironomus species. 
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Morphology (Table 6). Our specimens of C. staegeri are large sized plumosus-type larvae with a pale 
frontoclypeus. This latter feature clearly distinguishes them from C. ‘tigris ’ larvae that have a dark frontoclypeus. 
Outside of our study area, some C. staegeri larvae are reported to have a slightly darkened frontoclypeus (Martin 
2013). Thus other features such as the length of the lateral tubules and the outline of the anterior margin of 
ventromental plates could be examined to separate these two species. In our study area, C. staegeri could be 
distinguished from C. frommeri by the fact that the ventral and lateral tubules of the former were about half the 
length of those of the latter (Table 6). However, tubule length is not likely a reliable character to separate these 
species because in the study by Sublette and Sublette (1971) tubule lengths overlapped between these species. 
Sublette and Sublette (1971) suggested that C. frommeri and C. staegeri larvae could be separated by the structure 
of the anterior margin and apex of their paralabial plates as well as the shape of the teeth of the pecten 
epipharyngis. However, these features did not reliably separate these species in our study area. 

Cytology. The cytology of the single larva that we examined indicates that this species is C. staegeri since it 
has three chromosomes with a modified thummi arm combination of AB, CD, GEF (Wulker & Martin 1971; 
Kiknadze et al. 2004; Kiknadze et al. 2010). Thus it is distinct from C. ‘tigris ’, which possesses 2 chromosomes, 
and from C. frommeri, which has 4 chromosomes. 

Distribution and ecology (Table 7). C. staegeri has been found in a variety of lentic habitats from deep lakes 
to shallow pools (Wtilker et al. 1971) throughout Canada (British Columbia, Manitoba, Newfoundland, Northwest 
Territories, Ontario and Saskatchewan) and the United States (Alabama, California, Idaho, Illinois, Iowa, Kansas, 
Louisiana, Massachusetts, Michigan, Minnesota, Missouri, New Hampshire, New Jersey, New Mexico, New York, 
North Carolina, North Dakota, Pennsylvania, South Carolina, South Dakota, Tennessee, Washington and 
Wisconsin) (Sublette & Sublette 1971; Oliver et al. 1990; Martin et al. 2008; Martin 2013). In our study, C. 
staegeri was found in oligotrophic to eutrophic lakes and in a pond at depths ranging from 1-10 m, and at pH 
values ranging from 5.9-8.0. 

Taxonomic comment. Given the polymorphism of chromosomal inversions in populations of C. staegeri in 
Canada and the Unites States, Martin and Wulker (1971) speculated that C. staegeri might be in the process of 
splitting into three species based in part on their restriction to waters of different depths. Our DNA data do not 
support this idea since there is little variation in the coxl nucleotide sequences between C. staegeri that we 
collected from a pond, and over a range of depths in several lakes. In fact, the mean coxl intraspecific divergence 
among C. staegeri sequences is very low (0.04%). The different distributions of chromosomal inversions might 
therefore be due to populations with different inversion sequences adapting to different ecological niches, as has 
been suggested for species such as C. plumosus (Butler et al. 1999). 


Chironomus ( Chironomus ) frommeri Sublette and Sublette (1971) 

Material examined (Table SI): 7 larvae collected from an unnamed pond on a military base near Trois-Rivieres, 
Quebec. 

DNA (Figs. 8-9 and Tables 4-5, S3). In the coxl ID-tree, C. frommeri nucleotide sequences form a distinct 
clade. However, as mentioned above, coxl sequence divergences between C. frommeri and C. ‘tigris ’ (2-3%) as 
well as C. frommeri and C. staegeri (3-4%) are within the intraspecific sequence divergence range of Chironomus 
species assessed in this study (0-3%). Therefore, coxl sequence divergence values cannot be used to separate C. 
frommeri from C. ‘tigris ’ or C. staegeri. We were not successful in amplifying the gb2f gene for this species. 
Chosen restriction enzymes for the coxl PCR-RFLP analysis correctly separated C. frommeri from the other 
Chironomus species. 

Morphology (Table 6). Our specimens are large sized plumosus-type larvae with posterior ventral tubules 
longer than the anterior ones. Morphologically, larvae of C. frommeri from our study area can be distinguished 
from those of C. ‘tigris ’ and C. staegeri (however, see comments above on the morphology of these species). 

Cytology. The cytology of the larva analyzed permitted us to identify this species as C. frommeri. C. frommeri 
has four polytene chromosomes with the thummi arm combination of AB, CD, EF, G (Wulker & Martin 1971) as 
opposed to C. ‘tigris ' and C. staegeri that possess only 2 and 3 chromosomes, respectively. Arm G homologs of C. 
frommeri are closely paired, with a virtually terminal nucleolus, similar to the fused arm G of C. crassicaudatus 
(not collected in our study) (Wulker & Martin 1971). There is also another nucleolus proximal in arm B. 
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Distribution and ecology (Table 7). We were surprised to collect C. frommeri in eastern Canada because all 
previous collections of this species are from the western United States (California, Oregon, Utah and New Mexico) 
(Wtilker et al. 1971; Oliver et al. 1990). A possible explanation for this apparent anomaly is that this species was 
transported from the west to the east via military equipment since it was collected in a pond located in a military 
base. This species is known to occur in lakes, oxbows and permanent ponds (Martin 2013). 


Chironomus ( Chironomus ) cucini Webb (1969) 

Material examined (Table SI): 31 larvae collected in lakes in Quebec (Lake Bousquet, Lake Opasatica, Lake St. 
Joseph and Lake Vaudray) and in Ontario (Clearwater Lake). 

DNA (Figs. 8-9 and Tables 4-5, S3). We were able to obtain PCR products for both the coxl and gb2[S genes. 
In both of these ID-trees, C. cucini forms a distinct cluster. Consequently both genes can be used to accurately 
separate and identify C. cucini. The gb2p gene sequence of C. cucini has no intron, but contains three extra base 
pairs at the 3' end of the sequence. Chosen restriction enzymes for the coxl PCR-RFLP analysis correctly separated 
C. cucini from the other Chironomus species. 

Morphology (Table 6). Our specimens are large sized salinarius-type larvae lacking lateral and ventral tubules. 
In other study areas, larvae of C. cucini are reported to occasionally have small posterior ventral tubules (Wtilker & 
Butler 1983). Morphologically, larvae of C. cucini are similar to those of C. sp. NAIII with two minor differences. 
First, the structure of the pecten epipharyngis differs slightly between the species (Table 6). Second, the mean (± 
95% Cl) ratio of the lengths of antennal segments 1 to 2-5 (AR) of C. sp. NAIII larvae (1.77 ± 0.08) was 
significantly lower than that of C. cucini (2.04 ± 0.08). Note however that there was overlap in the ranges of the 
ARs between the two species. 

Cytology. The cytology of the 5 larvae analyzed permitted us to identity this species as C. cucini. C. cucini has 
four polytene chromosomes attached together by a chromocenter and with the thummi arm combination of AB, 
CD, EF, G (Martin 1979; Wtilker & Butler 1983). In most populations, a single nucleolus is located in arm G, 
although in some California populations there is a second nucleolus in arm B. 

Distribution and ecology (Table 7). C. cucini has been reported from across the Nearctic region (British 
Columbia, California, Indiana, Minnesota, New York and Ontario) (Oliver et al. 1990; Martin 2013). In our study 
area, C. cucini was found in the profundal zone of circum-neutral (pFl 6.2-7.5), oligotrophic to mesotrophic, lakes. 


Chironomus sp. NAII 

Material examined (Table SI): 4 larvae from Silver Lake (Ontario). 

DNA (Figs. 8-9 and Tables 4-5, S3). In the coxl ID-tree, C. sp. NAII nucleotide sequences clearly form a 
distinct clade. We were not successful in amplifying the gb2/3 gene for this species. Available DNA barcodes for 
Chironomus species did not cluster with those of C. sp. NAIL Coxl PCR-RFLP analysis correctly separated C. sp. 
NAII from the other Chironomus species. 

Morphology (Table 6). Our specimens are medium sized salinarius-type larvae. Morphologically, larvae of C. 
sp. NAII closely resemble those of C. cucini and C. sp. NAIII, but differ by the lobed dark spot in the middle of its 
frontoclypeus, by the partial coloration of its 3 rd mandibular tooth and by the type of mentum central trifid tooth. 

Cytology. This species has four polytene chromosomes with two nucleoli, one of which is in arm G. The 
cytology does not correspond to any known salinarius-type specimens from North America or Europe. The main 
difference between C. sp. NAII and the other known salinarius-type species is its lack of heterochromatic 
centromeres and its banding sequence in arm G. The cytological preparations were generally too poor to determine 
further details. This may be the larva of a previously described northern Chironomus species, for which the larva is 
currently unknown, or it may be a completely new species. 

Distribution and ecology (Table 7). This species was found in oligotrophic Silver Lake in 2007 at a depth of 4 
m. In 2010 and 2011, we sampled the lake again in an effort to collect additional C. sp. NAII larvae. Flowever, they 
were no longer present, perhaps because the pH of this lake had increased from 5.9 to 7.0 between 2007 and 2010. 
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Chironomus sp. NAIII (possibly C. decumbens (Malloch 1934)) 

Material examined (Table SI): 52 larvae collected in Lake D’Alembert (Quebec) and Hannah Lake, McFarlane 
Lake, Raft Lake and Ramsey Lake (Ontario). 

DNA (Figs. 8-9 and Tables 4-5, S3). We were able to obtain PCR products for both the coxl and gb2fi genes. 
C. sp. NAIII sequences form a distinct clade in both the coxl and gb2f ID-trees. Like C. cucini, the C. sp. NAIII 
gb2f sequence contains 3 extra base pairs at the 3' end, but unlike C. cucini, whose gb2[l sequence has no intron, 
the C. sp. NAIII gb2f sequence contains a type II intron. Additionally, unlike the other Chironomus species (except 
C. ‘tigris ’) the gb2/3 of C. sp. NAIII is three codons short immediately after the end of the 2 nd intron. Coxl PCR- 
RFLP analysis correctly separated C. sp. NAIII from the other Chironomus species. Sequences of C. sp. NAIII did 
not cluster with any of the available Chironomus species reference sequences. 

Morphology (Table 6). Our specimens are medium sized salinarius-type larvae that are difficult to distinguish 
from those of C. cucini (see comments on the separation of these species under C. cucini). The morphology of C. 
sp. NAIII (see Table 6) is similar to that of C. decumbens (see Martin 2013). 

Cytology. This species has three polytene chromosomes with heterochromatic centromeres. The arm 
combination is modified thummi-complex AB, CD, GEF. A nucleolus is located near the junction of arm G with 
arm E and a Balbiani ring is located towards the other end of arm G. Cytologically, this species fits the description 
of the North American cytospecies C. sp. 2x (Martin 2013) from Alaska, which is thought to be C. decumbens (Jim 
Sublette, personal communication). The only difference between our specimens and C. sp. 2x, is that our larvae 
possess a heavily heterochromatic centromere. This difference may or may not be significant since the presence of 
a heavily heterochromatic centromere can differ between populations and its detection can vary with the stain used. 
The voucher C. sp. 2x slide was stained with a brand of orcein that gave much paler staining. 

Unfortunately, no coxl or gb2f sequences of C. decumbens voucher specimens were available for comparison 
with sequences for our study larvae. Further investigation on a possible relationship between C. decumbens and C. 
sp. NAIII is clearly warranted. 

Distribution and ecology (Table 7). This species was found in oligotrophic to mesotrophic lakes at depths 
varying from 5-12 m and at pHs varying from 7.1-7.9. 


Chironomus sp. NAI (C. anthracinus- group) 

Material examined (Table SI): 9 larvae from Kasten Lake (Ontario). 

DNA (Figs. 8-9 and Tables 4-5, S3). We were able to obtain PCR products for both the coxl and gb2f genes. 
This species forms a distinct clade in the coxl ID-tree. However in the gb2/3 ID-tree, sequences of C. sp. NAI 
cluster with those of C. anthracinus. Chosen restriction enzymes for the coxl PCR-RFLP analysis correctly 
separated C. sp. NAI from the other Chironomus species. Available DNA barcodes did not allow us to identify C. 
sp. NAI (see discussion at the end of this section). 

Morphology (Table 6). Our specimens are large sized thummi-type larvae with the anterior ventral tubules 
slightly longer than the posterior ones. Morphologically, larvae of C. sp. NAI and C. anthracinus are 
indistinguishable. 

Cytology. C. sp. NAI has four short polytene chromosomes with the thummi arm combination of AB, CD, EF, 
G. The most common sequence in each chromosome arm is similar to that in C. anthracinus. 

Distribution and ecology (Table 7). This species was found in an oligotrophic, circum-neutral (pH 6.8), lake 
at a depth of 7.5 m. 

Taxonomic comment. C. sp. NAI larval morphology and cytology strongly resemble those of C. anthracinus. 
Additionally, gb2f sequences of C. sp. NAI and C. anthracinus cluster together in the ID-tree (sequence divergence 
varies from 1 to 5%). In fact, there are no consistent base differences between these two gb2fi sequences. Despite 
this lack of morphological, cytological or genetic difference, coxl sequences suggest that C. sp. NAI is a distinct 
species. The coxl sequence divergence between C. sp. NAI and C. anthracinus is relatively high (4-6%), and these 
sequences consistently differ by 22 bases (Table S4). For comparative purposes, other interspecific differences can 
be much lower, with C. staegeri and C. ‘tigris ’ differing by only 9 specific bases, C. frommeri and C. ‘tigris ’ by 13 
specific bases, and C. staegeri and C. frommeri by specific 19 bases. 
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In the coxl ID-tree, the reference sequence for Palearctic C. anthracinus larvae from Lake Esrom (Denmark) 
clusters with larvae from our lakes that have been cytologically identified as C. anthracinus, not with the adjacent 
cluster formed by our C. sp. NAI sequences. We hypothesize that larvae of C. sp. NAI might be C. rempelii 
Thienemannn (1941). Currently, there are conflicting opinions as to the status of C. rempelii. Based on adult 
morphology, Townes (1945) concluded that C. rempelii was a synonym of C. anthracinus Zetterstedt (1860). 
Shobanov et al. (1996) and Kiknadze et al. (2005) reached the same conclusion when they compared the 
chromosomes of these species. However, the heterochromatin on arm F and the sequences A3, C3 and F3 have so 
far been found only in samples from western Canada that include the type locality of C. rempelii (British Columbia, 
Alberta, Saskatchewan and Manitoba; Kiknadze et al. 2005). The large hetrochromatic block in the original C. 
rempelii population occurred in all males. However, no other populations were sexed, and the smaller blocks and 
the inversions are rare. Consequently, the absence of these in our material is not conclusive. Samples from the type 
locality of C. rempelii would be required to confirm our hypothesis. We amplified the coxl sequence of a 
Chironomus larva from British Columbia (Marion Lake) that was morphologically and cytologically 
indistinguishable from that of C. anthracinus (labelled as “C. {anthracinus-group.)" in Fig. 8) and found that its 
coxl sequence (Fig. 8) clusters with sequences of C. sp. NAI. Thus coxl sequences suggest that the currently 
recognized C. anthracinus in North America is not a single species, but a complex of at least two closely related 
species—the C. anthracinus-gmup. This may be the result of a recent speciation that occurred without 
hybridization, so that coxl has differentiated while the short and slower evolving gb2/3 sequence has not had time 
to accumulate significant changes. Further investigations of C. sp. NAI and the species status of C. rempelii are 
clearly warranted. 


Chironomus ( Chironomus) anthracinus Zetterstedt (1860) 

Material examined (Table SI): 27 larvae from Hannah Lake, Pine Lake, Raft Lake, Ramsey Lake and Silver Lake 
in Ontario and from Lake Arnoux and Lake Osisko in Quebec. 

DNA (Figs. 8-9 and Tables 4-5, S3). We were able to obtain PCR products for both the coxl and gblfi genes. 
C. anthracinus sequences form a distinct clade in the coxl ID-tree and cluster with the Palearctic C. anthracinus 
reference sequence, which confirms the identification of this species. In the gb2fl ID-tree, C. anthracinus 
sequences cluster with sequences of C. sp. NAI, the identity of which is uncertain (see discussion in the C. sp. NAI 
section). Chosen restriction enzymes for the coxl PCR-RFLP analysis correctly separated C. anthracinus larvae 
from the other Chironomus species. 

Morphology. Our specimens are large sized thummi-type larvae whose anterior ventral tubules are usually 
longer than their posterior ones. Morphologically, larvae of C. anthracinus and C. sp. NAI are identical. 

Cytology. The cytology of the 5 larvae analyzed was consistent with C. anthracinus, but identical to C. sp. 
NAI. C. anthracinus has 4 relatively short chromosomes with the thummi arm combination of AB, CD, EF, G 
(Kiknadze et al. 2005). The short chromosomes mean that banding patterns are often difficult to see clearly. There 
are two nucleoli: one on arm G and the other on arm F. 

Distribution and ecology (Table 7). C. anthracinus is widely distributed in the Holarctic region. In the 
Nearctic region, it occurs across Canada (Alberta, British Colombia, Manitoba, Ontario and Saskatchewan) and the 
United States (California, Indiana, Massachusetts, New Hampshire New York and Wisconsin) (Sasther 1975; 
Oliver et al. 1990; Saether 2012; Martin 2013). We found C. anthracinus at depths ranging from 2-12 m in highly 
acidic (pH 4.4) to circum-neutral (pH 8.5) lakes. The trophic status of these lakes was oligotrophic to mesotrophic, 
which is consistent with Saether’s (1975) suggestion that Nearctic C. anthracinus are more common in 
intermediate- to low-productivity lakes. In contrast, in the Palearctic region, C. anthracinus is known to be more 
frequently found in the profundal zone of moderately eutrophic lakes (Saether 1975). 


Chironomus ( Chironomus) entis Shobanov (1989) 

Material examined (Table SI): 8 larvae from Lake D’Alembert, Lake Dasserat, Lake Marlon. Lake Opasatica and 
Lake Pelletier (Quebec). 
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DNA (Figs. 8-9 and Tables 4-5, S3). We were able to obtain PCR products for both the coxl and the gb2p 
genes. In accordance with the findings of Makarevich et al. (2000), C. entis gh2p has a type II intron. 

In the coxl ID-tree, sequences of cytologically identified C. plumosus and C. entis larvae cluster together. In 
fact, some C. plumosus and C. entis sequences are identical. Similarities between the mitochondrial nucleotide 
sequences of these two species were also observed by Guryev and Blinov (2002), who found that trees based on the 
mitochondrial cytb gene did not group populations of C. entis and C. plumosus according to their species affiliation 
but rather according to their geographic occurrence. They attributed this phenomenon to mitochondrial gene flow 
that occurs when populations of sympatric sibling species produce fertile hybrids such that mitochondrial DNA 
appears in the progeny of the backcross. 

In contrast, in accordance with the findings of Guryev and Blinov (2002), our ID-tree based on nuclear gb2(3 
gene sequences successfully groups C. entis and C. plumosus according to species. The intraspecific variability of 
the partial gb2p nucleotide sequences of our C. entis (0%) and C. plumosus (2.3%) larvae are considerably lower 
that their interspecific variability (15-16%). Thus C. entis and C. plumosus can be distinguished using the gb2p 
gene, but not the coxl gene. Lastly, since PCR-RFLP analysis was performed using the coxl gene, this result could 
not be used to separate C. entis and C. plumosus. 

Morphology (Table 6). Our specimens are very large semireductus-type larvae. Morphologically, C. entis and 
C. plumosus are almost indistinguishable. Kiknadze et al. (1991) described the outer hooks on the anterior margin 
of the ventromental plates as being shorter and blunter in C. plumosus than in C. entis in Palearctic populations; 
however, we did not observe such differences in our specimens. C. entis are semireductus-type larvae, whereas 
those of C. plumosus vary from being semireductus-type to plumosus-type. In fact, some C. plumosus larvae from 
our study lakes have ventral tubules that are intermediate between the plumosus-type and semireductus-type. 
Within a given lake, the ventral tubules of C. plumosus were always longer than those of C. entis, which allowed 
the larvae of these species to be separated. However, between lakes, there was considerable overlap in their lengths. 
Consequently, C. plumosus and C. entis cannot be distinguished based solely on morphology such that cytological 
and genetic techniques are needed to unambiguously separate them. 

Cytology. All larvae were analyzed cytologically to confirm the identification of this species and to verily the 
accuracy of the DNA results. C. entis has four relatively short chromosomes with the thummi arm combination of 
AB, CD, EF, G (Kiknadze et al. 2000a; Kiknadze et al. 2000b; Gunderina et al. 2009), with only a single nucleolus 
in arm G. It shares two rare sequences with Nearctic populations of C. plumosus (Kiknadze et al. 2000a). 

Distribution and ecology (Table 7). C. entis has previously been reported from lakes in Canada (British 
Columbia, Manitoba, Ontario and Saskatchewan) and in the United States (Colorado, Indiana, Minnesota, North 
Dakota, Oklahoma, South Dakota and Wisconsin) (Kiknadze et al. 2000a; Martin 2013). We found C. entis in 
mesotrophic to eutrophic, circum-neutral (pH 7.1-8.3) lakes at depths ranging from 1-9 m. 


Chironomus ( Chironomus ) plumosus Linnaeus (1758) 

Material examined (Table SI): 33 larvae from lakes in Quebec (Lake D’Alembert, Lake Duprat, Lake Fortune, 
Lake Kinojevis, Lake Marlon, Lake Osisko, Lake Pelletier, Lake Rouyn and Lake Saint Augustin) as well as Kelly 
Lake in Ontario. 

DNA (Figs. 8-9 and Tables 4-5, S3). We were able to obtain PCR products for both the coxl and gb2p genes. 
In accordance with the findings of Makarevich et al. (2000), C. plumosus gb2p\\ , && a type II intron. C. plumosus 
cannot be distinguished from C. entis on the basis of coxl sequences. However, these species can be separated 
using gb2p sequences (see comments in the C. entis DNA section). 

Morphology (Table 6). Our specimens are very large semireductus (commonly referred to as “semi-reductus” 
type in reference to C. plumosus )—to plumosus-type larvae (see comment in section on C. entis morphology). 
Morphologically, larvae of C. entis and C. plumosus are indistinguishable (see comment in C. entis morphology 
section). 

Cytology. All larvae were analyzed cytologically to confirm the identification of this species and to verily the 
accuracy of the DNA results. C. plumosus has four relatively short chromosomes with the thummi arm 
combination of AB, CD, EF, G (Butler et al. 1999) with only a single nucleolus in arm G. 

Distribution and ecology (Table 7). This species was previously known from lakes in Canada (British 
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Columbia, Manitoba, Ontario and Saskatchewan) and the United States (Alabama, California, Colorado, Indiana, 
Kentucky, Massachusetts, Minnesota, New Mexico, North Dakota, Oklahoma, South Dakota and Wisconsin) 
(Butler et al. 1999; Martin 2013) at depths up to 23m (Martin 2013). We found C. plumosus at depths ranging from 
1-8 m, and in oligotrophic to eutrophic lakes ranging in pH from 6.8-8.5. 


Chironomus ( Chironomus ) maturus Johannsen (1908) 

Material examined (Table Si): 10 larvae from Lake Bedard (Quebec). 

DNA (Figs. 8-9 and Tables 4-5, S3). We were able to obtain PCR products for both the coxl and the gb2ft 
genes. In both ID-trees, C. maturus sequences form a distinct clade. In the coxl ID-tree, sequences of the larvae we 
collected cluster with the reference sequence of C. maturus, thereby confirming the identification of this species. 
Coxl PCR-RFLP analysis accurately separated C. maturus larvae from the other Chironomus species. 

Morphology (Table 6). Our specimens are large sized plumosus-type larvae with long ventral and lateral 
tubules. The frontoclypeus is pale; however, some C. maturus from other regions are reported to have a dark 
frontoclypeus (Martin 2013). Larvae of C. maturus could be distinguished from the other collected plumosus-type 
larvae through a combination of morphological features (see morphological key at the end of this section). 

Cytology. C. maturus larvae were not identified through cytology because their chromosomes were not in 
good enough condition. C. maturus is known to possess four polytene chromosomes with a maturus arm 
combination of AF, BE, CD, G. The cytology of C. maturus has been described by Wiilker and Martin (1974) and 
Kiknadze et al. (2004). 

Description and ecology. This species has been recorded previously from shallow pools and polluted water 
bodies (Martin 2013) in central Canada (Manitoba and Ontario) and the United States (Alaska, California, Dakota, 
Indiana, Louisiana, New Mexico, New York, South Dakota and Wisconsin) (Oliver et al. 1990; Martin 2013). We 
collected C. maturus in a mesotrophic and circum-neutral (pH = 7.3) lake. 


Chironomus ( Chironomus) tlecorus -group sp. 2 Butler et al. (1995) 

Material examined (Table SI): 17 larvae from Lake Adeline, Lake Dufault, Lake Duprat, Lake Fortune, Lake 
Opasatica and the St. Charles River in Quebec as well as Silver Lake in Ontario. 

DNA (Figs. 8-9 and Tables 4-5, S3). In the coxl ID-tree, C. decorus- group sp. 2 nucleotide sequences form a 
distinct clade. Sequences of larvae we collected cluster with the reference sequence of C. decorus- group sp. 2 from 
GenBank, therefore confirming the identity of this species. The coxl interspecific sequence divergence between C. 
decorus- group sp. 2 and reference sequence C. quinnitukqut (3%) (data not shown) is within the intraspecific 
sequence divergence range (0-3%) of Chironomus species assessed in this study. Therefore, sequence divergence 
cannot be used to separate these species. Coxl PCR-RFLP analysis correctly separates C. decorus- group sp. 2 from 
the other Chironomus species. We were not successful in amplifying the gb2p gene of this species. 

Morphology (Table 6). Our specimens are mostly medium sized bathophilus- or fluviatilis-type larvae (Table 
SI). However, whereas Quebec specimens had no lateral tubules, those from Silver Lake, Ontario, had small lateral 
tubules (melanotus-type). The 3 rd inner tooth of the mandibles is pale and fused to the lower margin, however, in 
lakes from other regions, they are reported to be partially darkened (Martin 2013). From the morphological 
description of C. quinnitukqut in Martin (2013), larvae of C. decorus- group sp. 2 that do not possess lateral tubules 
could not be distinguished from those of C. quinnitukqut. Furthermore, larvae of C. decorus- group sp. 2 cannot be 
distinguished from those of C. bifurcatus. 

Cytology. The cytology of the two larvae analyzed clearly indicates that they belong to a species currently 
referred to as C. decorus- group sp. 2. C. decorus-gmwp sp. 2 has four polytene chromosomes with the thummi arm 
combination of AB, CD, EF, G (Butler et al. 1995). Typical of members of the C. decorus- group, it has only a 
single nucleolus, which is virtually terminal in arm G. 

Description and ecology. This species has been collected previously in Canada (Saskatchewan) and the 
United States (Massachusetts, Minnesota, Mississippi, New Mexico, North Dakota, Vermont and Wisconsin) 
(Butler et al. 1995; Martin 2013). C. decorus- group sp. 2 larvae are reported from depths greater than 10 m (Martin 
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2013); however, in our study C. decorus -group sp. 2 larvae were found at depths of 1-5 m. C. decorus -group sp. 2 
was collected in oligotrophic to mesotrophic lakes of circum-neutral pH (7.2-7.6). However, it was also found in 
the St. Charles River where sediments have been contaminated by untreated municipal waste waters. 


Chironomus ( Chironomus ) harpi Sublette (in Wiilker et al. 1991) 

Material examined (Table SI): 6 larvae from Lake Arnoux (Quebec). 

DNA (Figs. 8-9 and Tables 4-5, S3). We were able to obtain PCR products for both the coxl and gb2j3 genes. 
C. harpi sequences form distinct clades in both the coxl and gb2f ID-trees. Coxl sequences of collected larvae 
cluster with the reference sequence of C. harpi, thus confirming the identification of this species. Coxl PCR-RFLP 
profiles were not obtained for C. harpi. However, through simulation digests, we were able to demonstrate that the 
chosen restriction enzymes (SspI, Hinfl , Rsal and Taql) would have accurately separated C. harpi from the other 
Chironomus species. 

Morphology (Table 6). Our specimens fit the morphological description of C. harpi (Martin 2013). They are 
medium sized plumosus-type larvae with posterior ventral tubules longer than the anterior ones. Morphologically, 
larvae of C. harpi strongly resemble those of the other plumosus-type larvae collected in our study by having a pale 
frontoclypeus (C. frommeri, C. staegeri and C. maturus). C. harpi larvae can be distinguished from larvae of these 
other species principally by the type of the middle trifid tooth on the mentum (see Table 6). 

Cytology. We did not identify C. harpi through cytology. They are reported to possess four polytene 
chromosomes with the thummi arm combination of AB, CD, EF, G, with a large nucleolus near the centromere of 
arm D, and a second nucleolus sometimes developed medially in arm G. (Wiilker et al. 1991). 

Distribution and ecology (Table 7). In our study, we found C. harpi at depths of 1-4 m in a lake that has been 
heavily impacted by acid mine drainage (pH 2.7-3.8). This further supports the identification of this species, since 
C. harpi had been reported previously only from acidic pools (Arkansas, Illinois, New York and South Dakota) 
(Martin 2013). 

Chironomus ( Chironomus) bifurcatus Wuelker, Martin, Kiknadze, Sublette and Michiels (2009) 

Material examined (Table SI): 50 larvae collected in Quebec (Lake Adeline, Lake Arnoux, Lake D’Alembert, 
Lake Duprat, Lake Kinojevis, Lake Opasatica and Lake St. Joseph) and in Ontario (McFarlane Lake and Silver 
Lake). 

DNA (Figs. 8-9 and Tables 4-5, S3). We were able to obtain PCR products for both the coxl and gb2f genes. 
In the coxl ID-tree, sequences of collected larvae cluster with the reference sequences of C. bifurcatus, thereby 
confirming the identification of this species. However, in the coxl ID-tree, collected and reference C. bifurcatus 
nucleotide sequences form two distinct clades (labelled in Fig. 8 as group 1 and group 2) that differed by 7 specific 
bases (Table S5). There were also consistent cytological differences between these two groups (see cytology 
section). However, the sequence divergence between coxl -sequences from these two groups is low (2%) and, in the 
gb2[ J > ID-tree, C. bifurcatus sequences from both groups cluster together. 

Three different profiles were created through PCR-RFLP analysis of the coxl gene. One profile included all 
group 1 larvae and the other two profiles (obtained using the Taql enzyme) included all group 2 larvae. Simulation 
digests demonstrate that the SspI, Hinfl, Rsal and Taql restriction enzymes would not have separated C. bifurcatus 
(gr. 2) larvae from C. dilutus because in some cases their restriction fragment lengths are so similar that it would be 
difficult to differentiate them on an agarose gel. However, simulation digests demonstrate that both species should 
be separable using the restriction enzyme Alul. 

Morphology (Table 6). Our specimens are medium sized bathophilus-type larvae with the anterior ventral 
tubules slightly longer than the posterior ones. All of our larvae lacked lateral tubules. In contrast, larvae from 
some populations of C. bifurcatus are reported to have either very small (180 pm) lateral tubules (about 180 pm; 
Wuelker et al. 2009) (melanotus-type) or be fluviatilis-type (Martin 2013). The frontoclypeus of our C. bifurcatus 
larvae is pale, but is reported to be slightly darkened in larvae from some other regions (Wuelker et al. 2009). There 
were no morphological differences between larvae belonging to groups 1 and 2. Morphologically, larvae of C. 
bifurcatus cannot be distinguished from those of C. decorus -group sp. 2. 
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Cytology. The cytology of the 5 specimens analyzed permitted us to identity them as C. bifurcatus. Larvae of 
this species have four polytene chromosomes with the thummi arm combination of AB, CD, EF, G, as well as a 
single, virtually terminal, nucleolus in arm G (Wuelker et al. 2009). Flowever, we noted some cytological 
differences among our larvae and among reference specimens that allowed all of the specimens to be separated into 
two groups that corresponded to those mentioned above for the coxl gene. Thus larvae from group 1 have 
cytological sequences B1 and FI with no median Balbiani ring in the middle of arm G whereas those from group 2 
have cytological sequences B2 and F2, and a Balbiani ring in the middle of arm G. 

Distribution and ecology (Table 7). This species has been collected previously in southern Canada (Manitoba, 
Ontario and Quebec) and the northern United States (Massachusetts, Michigan, Minnesota and Wisconsin) 
(Wuelker et al. 2009; Martin 2013). We collected larvae of C. bifurcatus gr. 1 only in Quebec lakes (Lake Arnoux, 
Lake D’Alembert, Lake Duprat, Lake Opasatica and Lake St. Joseph), whereas those in gr. 2 were found in lakes in 
both Quebec (Lake Kinojevis and Lake St. Joseph) and Ontario (McFarlane Lake and Silver Lake). Larvae 
belonging to the two genetic groups differed somewhat in their distribution with respect to lake water pH, trophic 
status and water depth. Specifically, larvae of C. bifurcatus (gr. 1) were collected in acidic to circum-neutral (pH 
2.7-7.6) and oligotrophic to eutrophic lakes at depths ranging from 1.5-9 m. In contrast, larvae of C. bifurcatus (gr. 
2) were collected in circum-neutral (pH 7.0-7.8) and oligotrophic to mesotrophic lakes over a wider range of 
depths (4-24m). We note that larvae of C. bifurcatus (gr. 2) were always found along with those of C. staegeri and 
C. ‘tigris ’ whereas this was not the case for gr. 1 larvae. 

Taxonomic comments. Although coxl sequences and cytological differences separate C. bifurcatus into two 
groups, we suggest that it is premature to recognize two species until further work has been completed. 


Chironomus ( Chironomus ) dilutus Shobanov, Kiknadze and Butler (1999) 

Material examined (Table SI): 9 larvae from Kelly Lake (Ontario). 

DNA (Figs. 8-9 and Tables 4-5, S3). We were able to obtain PCR products for both coxl and gh2fi genes. In 
accordance with the findings of Makarevich et al. (2000), C. dilutus gb2/3 has no intron. In the coxl ID-tree (Fig. 
8), reference sequences of C. pallidivittatus (sensu Beermann 1955) cluster with reference sequences of C. dilutus. 
Similarities between the mitochondrial nucleotide sequences of these two species were also observed by Martin et 
al. (2002), who showed that trees based on mitochondrial sequences clustered populations of C. dilutus and C. 
pallidivittatus according to their geographic distribution, whereas those based on nuclear sequences clustered 
populations according to their species affiliation. The inability of the mitochondrial coxl gene to separate C. 
dilutus and C. pallidivittatus is likely due to mitochondrial gene flow (Martin et al. 2002). To confirm the 
identification of larvae whose coxl sequences clustered with reference C. dilutus and C. pallidivittatus sequences, 
we amplified and sequenced the partial gb2f gene. Partial gb2p gene sequences of all larvae were identical to the 
C. dilutus reference sequence (Fig. 9), which suggests that our larvae are not C. pallidivittatus but more likely C. 
dilutus. Coxl PCR-RFLP profiles were not obtained for C. dilutus. However, simulation digests demonstrated that 
the SspI, 11 infI. Rsal and TaqI restriction enzymes would not have separated C. dilutus from C. bifurcatus and C. nr. 
atroviridis (sp. 2i) because their restriction fragment lengths are so similar in some cases that it would be difficult 
to differentiate them on an agarose gel. However, simulation digests demonstrate that C. dilutus could be 
distinguished from C. bifurcatus and C. nr. atroviridis (sp. 2i) using the restriction enzyme A hiI. 

Morphology (Table 6). Our specimens are large sized plumosus-type larvae with their posterior ventral 
tubules usually longer than their anterior ventral tubules. Although the teeth of the mentum and mandibles of C. 
dilutus are reported to be rounded (Martin 2013), in our larvae, they varied from being rounded to sharp, which 
suggests that this feature cannot be used to identify larvae of C. dilutus in our study area. Morphologically, larvae 
of C. dilutus cannot be distinguished from those of C. pallidivittatus (not collected in our study; see Martin 2013). 

Cytology. C. dilutus larvae have four polytene chromosomes with the camptochironomus arm combination of 
AB, DE, CF, G. The cytology of C. dilutus has been described by several authors (see Martin 2013). 

Distribution and ecology (Table 7). This species has been found previously in numerous localities across 
Canada (Alberta, British Columbia, Manitoba, Ontario and Saskatchewan) and the northern United States (Iowa, 
Massachusetts, Michigan, Minnesota, New York, North Dakota, South Dakota, Utah, Wisconsin and Wyoming) 
(Shobanov et al. 1999; Martin 2013). C. dilutus is known to thrive in organically-enriched eutrophic water bodies 
(Townes 1945, referred to as Tendipes (Tendipes) teutons (Fabricius)). We collected large numbers of C. dilutus in 
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eutrophic Kelly Lake (at 5 m, pH 7.5) where sediments have been highly contaminated by untreated sewage from 
the city of Sudbury (1880s—1972) along with discharges from mining, milling and smelting operations (1880s- 
present) (City of Greater Sudbury 2013). Most of the larvae we collected had deformed mouthparts, likely due to 
the contaminants to which they were exposed (Hare & Carter 1976). 


Chironomus ( Chaetolabis ) ochreatus Townes (1945) 

Material examined (Table SI): 2 larvae from Lake Opasatica (Quebec). 

DNA (Figs. 8-9 and Tables 4-5, S3). In the coxl ID-tree, C. ochreatus nucleotide sequences clearly form a 
distinct clade. Sequences of collected species cluster with the reference sequence of C. ochreatus thereby 
confirming the identification of this species. We were not successful in amplifying the gb2fi gene for this species 
and coxl PCR-RFLP profiles were not obtained for this species. However, simulation digests using the restriction 
endonucleases Sspl, Hinfi, Rsal and Taq I demonstrate that coxl PCR-RFLP analysis would have separated C. 
ochreatus from the other Chironomus species. 

Morphology (Table 6). Our specimens are medium sized thummi-type larvae. The 3 rd inner teeth of the 
mandibles are partially darkened and fused to the lower margin (type B; Fig. 6). In contrast, in specimens studied 
by Martin (2013), the 3 rd inner tooth was pale and well separated from the lower margin. 

Cytology. The cytology of the two larvae analyzed clearly indicates that this species is C. ochreatus. C. 
ochreatus has three polytene chromosomes that are thought to have a modified thummi arm combination of AB, 
CD, GEF (Martin 2012). Arm G is generally unpaired with a nucleolus near the junction with arm E. 

Distribution and ecology (Table 7). This species was previously known from the eastern United States 
(Arkansas, Maine, Massachusetts, Michigan, New Jersey, New York, Rhode Island, South Carolina, Virginia and 
Wisconsin) (Townes 1945; Oliver et al. 1990; Martin 2013). In our study, we collected C. ochreatus in a single 
mesotrophic Quebec lake at a depth of 2 m (pH 7.7). 


Chironomus {Chaetolabis) nr. atroviridis (sp. 2i) Martin (2013) 

C. atroviridis Townes has been found to comprise two species in North America, one with four polytene 
chromosomes (sp. 2i of Martin 2013), and the other with only three polytene chromosomes (sp. 2h of Martin 2013). 
Only the former species occurred in our samples. 

Material examined (Table SI): 4 larvae from Lake Marlon (Quebec). 

DNA (Figs. 8-9 and Tables 4-5, S3). In the coxl ID-tree, C. nr. atroviridis (sp. 2i) nucleotide sequences form 
a distinct clade. Sequences of collected larvae cluster with the reference sequence of C. nr. atroviridis (sp. 2i), 
which confirms the identification of this species. We were not successful in amplifying the gb2fi gene of this 
species. Coxl PCR-RFLP profiles were not obtained for this species. However, simulation digests demonstrate that 
Sspl, Hinfi, Rsal and Taql could not be used to separate C. nr. atroviridis (sp. 2i) from all C. dilutus specimens 
because their restriction fragment lengths are so similar that it would be difficult to differentiate them on an agarose 
gel. However, simulation digests demonstrate that both species could be separated using the restriction enzyme 
Alul. 

Morphology (Table 6). Our specimens are large sized thummi-type larvae. They are morphologically very 
similar to those of C. ochreatus, in that both species have the teeth of the pecten epipharyngis flattened, however, 
the gula of C. nr. atroviridis (sp. 2i) is darkened posteriorly whereas that of C. ochreatus is pale or only slightly 
darkened. 

Cytology. The cytology of the 2 larvae analyzed clearly indicates that this species is C. nr. atroviridis (sp. 2i). 
C. nr. atroviridis (sp. 2i) has four polytene chromosomes with some indication of a thummi arm combination 
(Martin 2013). Arm G is generally unpaired with a virtually terminal nucleolus. There are no nucleoli in the other 
chromosomes. 

Distribution and ecology (Table 7). This species has been collected previously in Manitoba and Ontario in 
shallow water near macrophytes (Martin 2013). Likewise, we collected C. nr. atroviridis (sp. 2i) in the vegetated 
littoral zone (1 m depth) of a mesotrophic to eutrophic lake (pH 7.4). 

Taxonomic comment. Wiederholm (1979) considered C. ochreatus to be a synonym of C. atroviridis, but 
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mentioned that further study was needed. Because Wiederholm (1979) was not aware that there were two forms of 
C. atroviridis (2i and 2h; Martin 2013), we do not know which form would correspond to the material he was 
comparing. In any case, our results clearly indicate that C. nr. atroviridis (sp. 2i) is distinct from C. ochreatus. First, 
the cox 1 sequences of C. nr. atroviridis (sp. 2i) and of C. ochreatus are strikingly different such that their average 
interspecific divergence (12%; Table 5) is much higher than their intraspecific divergences (<3%; Table S3). 
Second, C. nr. atroviridis (sp. 2i) has 4 chromosomes whereas C. ochreatus has only 3. C. ochreatus also differs 
from the three chromosome forms of C. nr. atroviridis (sp. 2h) by the position of the nucleolus, which is 
subterminal in the latter species. 


Morphological key to larvae of the Chironomus species collected in our study 

The following key for identifying fourth-instar larvae of the Chironomus species collected in our study is based on 
the morphology of the tubules on the 10 th and 11 th body segments and features of the head capsule. Illustrations of 
these features are given in Figures 1—7. 

Although the majority of our study species can be separated using the following key, we acknowledge that, as 
with most morphological classifications, the characters used are likely to show some variability due to genetic or 
environmental factors. In this key, larvae are first separated based on the absence/presence, length and shape of 
tubules, which was effective in separating larvae of all of our collected species with the exception of C. decorus- 
group sp. 2 since some larvae of this species had lateral tubules whereas others did not. The ventral tubules of C. 
decorus -group sp. 2 larvae also varied from being straight (bathophilus-type) to slightly curved (fluviatilis-type). 
Our C. bifurcatus larvae were all bathophilus-type, but in other geographical regions have been ascribed to several 
larval types (bathophilus, fluviatilis or melanotus; Martin 2013) and these differences are reported to be related to 
the depth or the type of substrate on which larvae occur (Martin 2013). At an extreme, the presence or absence of 
lateral tubules among C. bifurcatus larvae is reported to vary among larvae hatched from the single egg mass from 
which the type was reared. (J. Martin, unpublished). 

Note that this key is based on morphological features of the Chironomus species that we collected, such that 
other species in the Nearctic could fit these descriptions (see notes at the end of the key). 


1 

2 


3 

4 

5 

6 
7 


9 


10 


11 01 segment without ventral tubules; 10 th segment without lateral tubules (salinarius-type larvae). 2 

11 th segment with one or two pairs of ventral tubules.3 

Frontoclypeus with a dark longitudinal stripe and a lobed dark spot in the middle; central trifid tooth of mentum with outer 
teeth almost completely separated from middle tooth (type C); 3 rd inner tooth of mandibles partially darkened and fused to 

lower margin (type B) .C. sp. NAII 1 

Frontoclypeus pale or slightly darkened with a lobed dark spot in anterior portion; central trifid tooth of mentum with outer 
teeth partially separated from middle tooth (type B); 3 rd inner tooth of mandibles pale and fused to lower margin (type A) ... 

. C. cucini (in part) or C. sp. NAIII land2 

11 th segment with only one pair of short ventral tubules (located posteriorly); 10 111 segment without lateral tubules (halophilus- 


type) . C. cucini (in part) 

11 01 segment with two pairs of ventral tubules. 4 

10 th segment without lateral tubules.5 

lO 01 segment with a pair of lateral tubules.8 

Ventral tubules straight (bathophilus-type) or slightly curved (fluviatilis-type) . 

. C. bifurcatus (in part) or C. decorus- group sp. 2 (in part) 3 

Anterior ventral tubules with an elbow; posterior ventral tubules coiled (thummi-type). 6 

Pecten epipharyngis teeth flattened (type C or D). 7 

Pecten epipharyngis teeth elongated (type A or B) . C. anthracinus or C. sp. NAI 

Gular region almost completely pale or at most slightly darkened. C. ochreatus 

Gular region darkened posteriorly . C. nr. atroviridis (2i) 

Ventral tubules straight or slightly curved.9 

Anterior pair of ventral tubules with an elbow; posterior pair of ventral tubules coiled (plumosus-type larvae).10 

Ventral tubules equal to or greater than the width of the 11 th segment (melanotus type); 3 rd inner tooth of mandible pale. 

. C. decorus- group sp. 2 (in part) or C. bifurcatus (in part) 


Ventral tubules less than the width of the ll lh segment (semireductus-type larvae); 3 rd inner tooth of mandible dark. 

. C. ends or C. plumosus (in part) 


3 rd inner tooth of mandible partially dark to dark.11 

3 rd inner tooth of mandible pale. 12 
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11 

12 

13 

14 


15 


Frontoclypeus pale and gula strongly to completely darkened. C. plumosus (in part) 

Frontoclypeus dark and gula slightly to posteriorly darkened. C. dilutus 4 

Anterior margin of ventromental plates crenulated. C. frommeri or C. staegeri 5 

Anterior margin of ventromental plates smooth to relatively smooth.13 

Frontoclypeus pale.14 

Frontoclypeus dark. 15 

3 rd inner tooth of mandible fused to lower margin, central trifid tooth of mentum with outer teeth only partially separated from 
middle tooth (type B) and 4 th lateral teeth reduced to the height of the S” 1 lateral teeth (type II); found in highly acidic waters. . 

. C. harpi 

3 rd inner tooth of mandible separated from lower margin, central trifid tooth of mentum with outer teeth distinctly separated 

from middle tooth (type D) and 4 th lateral teeth only slightly reduced (type I). C. maturus (in part) 

Gula pale to slightly darkened and 4 th lateral teeth of mentum only slightly reduced (type I). C. maturus (in part) 

Gula strongly to completely darkened and 4 tt lateral teeth of mentum about the same height as the 5 th lateral teeth (type II) . . . 
. C. ‘tigris ’ 


Notes: 

It is likely that C. atritibia (Malloch 1934) would also key here, as it is reported to have a salinarius-type larva (Wiilker & But¬ 
ler 1983). Although C. atritibia is thought to have a more northerly distribution, we cannot rule out the possibility that it corre¬ 
sponds to either C. sp. NAII or C. sp. NAII1. 

In our specimens, the mean AR of C, cucini (2.08) was significantly greater than that of C, sp. NA1II (1.77). Although there 
was some overlap in the range of ARs between the two species, 4 of 5 C. cucini larvae had an AR >2.0, whereas 10 of 11 spec¬ 
imens of C. sp. NAIII had an AR <1.95.The more southern C. major (Wiilker & Butler 1983) would also key here, but is much 
larger (30-55 mm; Epler 2001). 

Other Nearctic species of the C. decorus- group (Wuelker 2010; Saether 2012) are also likely to key out here (ex. C. quinni- 
tukqut) 

4 C. pallidivittatus would also key out here. 

5 C. crassicaudatus (Malloch 1915) would also key out here. 


Approaches used to delimit Chironomus species 

No single approach (morphology, cytology, genetics) was adequate for delimiting and identifying larvae of all of 
the Chironomus species that we collected. 

Thus larval morphology alone could not be used to separate five pairs of Chironomus species that we collected, 
i.e., C. cucini and C. sp. NAIII, C. bifurcatus and C. decorus- group sp. 2, C. anthracinus and C. sp. NAI, C. 
staegeri and C. frommeri, as well as C. entis and C. plumosus. Furthermore, some Nearctic species that we did not 
collect are reported to be morphologically identical to our study species (see notes at the end of the morphological 
key). In addition, for some of our species, specimens from other regions are reported to differ morphologically 
from those that we collected, and it is known that some larval characters are affected by wear, environmental 
conditions and genetic variation (Martin 2013). This said, larval morphology was undeniably important when used 
in combination with other methods of species delimitation. 

Differences in larval cytology, based on the structure of salivary-gland polytene chromosomes, allowed the 
definitive identification of many of the species that we collected. Indeed, the cytology of most North American 
Chironomus species has been described (Martin 2013). However, one must bear in mind that cytology is faced with 
the same challenges as the other identification methods; that is, it is not always possible to determine whether or 
not differences in chromosome banding patterns and other structures are attributable to species differences or to 
regional or individual differences within a given species (Martin 2011). Indeed, only a handful of taxonomists 
worldwide have the necessary expertise to identify Chironomus species through cytology, which is a major 
drawback for non-cytological experts wishing to identify Chironomus species. 

Genetic techniques, namely PCR-RFLP analysis and DNA barcoding of the coxl gene, successfully separated 
and identified most of the Chironomus species that we collected. We present the first coxl sequences for many of 
the known Nearctic Chironomus species (C. acidophilus, C. anthracinus, C. bifurcatus, C. cucini, C. frommeri, C. 
harpi, C. nr. atroviridis (sp. 2i), C. ochreatus, C. plumosus, C. quinnitukqut, C. sp. g, C. sp. h and C. 'tigris ’). 
However, DNA barcoding failed to distinguish between two species pairs (C. entis and C. plumosus', C. dilutus and 
C. pallidivittatus) because each pair has identical coxl nucleotides sequences. Such sequence similarities are likely 
the result of mitochondrial gene flow, and have been found in a number of closely related species groups around the 
world (e.g. Martin 2011). 


432 ■ Zootaxa 3741 (4) © 2013 Magnolia Press 


PROULX ETAL. 













Using the gb2f gene, we successfully separated C. entis from C. plumosus and C. dilutus from C. 
pallidivittatus and confirmed the species statuses of C. cucini, C. bifurcatus, C. harpi, C. maturus, C. sp. NAIII and 
C. ‘tigris ’. We present the first published gb2fi sequences for C. anthracinus, C. bifurcatus, C. calligraphus, C. 
cucini, C. harpi, C. maturus, C. sp. u and C. ‘ tigris ’. One downside of using the gb2f gene was that we were not 
able to obtain PCR products for all species. In fact, no primer combination was able to amplify all Chironomus 
species, which is a limitation when using this gene for DNA barcoding. For these species we sometimes obtained 
two PCR products or sequences belonging to the globin 7A or 9 genes. Our results confirm those of Hankeln el al. 
(1997) who found that this gene is highly variable and that the only conserved regions are also conserved in the 
gb7A and gb9 genes. This is a major drawback for using the gb2[j gene for identifying Chironomus species. Thus, 
the use of another nuclear gene might be more appropriate for Chironomus species identifications. Other studies 
have used the nuclear internal transcribed spacer (ITS) region to separate C. plumosus from C. entis as well as other 
Chironomus species (Gunderina & Katokhin 2011; Martin 2011; Gunderina 2012). The nuclear 
carbamoylphosphate synthetase (CAD) region has also been successfully used to identify chironomid species 
(Carew et al. 2011). Another disadvantage of using the gb2f gene is that nuclear genes evolve more slowly than do 
mitochondrial genes. Thus, relatively recent speciation might not always be detected when using the gb2f gene. 
This might be the reason why sequences of C. sp. NAI and C. anthracinus, as well as those of C. bifurcatus (gr. 1) 
and C. bifurcatus (gr. 2), differ for the coxl gene but are identical for the gb2fi gene. Thus, incorporation of both 
mitochondrial and nuclear genes, whose modes of inheritance and mutation rate differ, clearly provides better 
resolution for Chironomus species identification. 

Several studies have advocated the use of sequence divergence thresholds to separate species (ex. Hebert et al. 
2004b). However, our results demonstrate that sequence divergence thresholds cannot be used to separate all 
Chironomus species. Thus we recorded overlap between intra- and inter-specific sequence divergences for both of 
the genes that we studied (Table 5); similar overlaps have been reported for other Chironomus species (Martin 
2011), as well as for species of other chironomid genera (Carew et al. 2005; Ekrem et al. 2007) and other types of 
dipterans (Meier et al. 2006). In our study, Chironomus coxl intraspecific sequence divergences were < 3% (Table 
S3). For most species, coxl interspecific sequence divergences ranged from 9% to 20%, but between some of our 
study species the interspecific divergences ranged from 1 to 4% (Table 5). This overlap was due either to some of 
our study species sharing identical sequences (i.e., C. entis/C. plumosus and C. dilutus/C. pallidivittatus ) or to 
interspecific divergences being so low that they fell within the intraspecific range for the genus (i.e., C. staegerilC. 
‘tigris 7C. frommeri (Table 5) and C. decorus- group sp. 2/C. quinnitukqut (data not shown)). With respect to the 
gb2f> gene, Chironomus species intraspecific sequence divergences were < 2%. The interspecific divergences 
between most species ranged from 5 to 46% (Table 5), but sequence divergences between C. sp. NAI and C. 
anthracinus ranged from 1 to 5% (Table 5), which is within the intraspecific range for species. In light of our 
results, the calculated intraspecific sequence divergences of 3% for the coxl gene and 2% for the gb2(3 gene can be 
used as a guide to help sort Chironomus species, but should not be used in isolation. DeSalle et al. (2005) have 
suggested that, rather than looking at sequence divergences, specific base differences that characterize related 
species should be sought. We used this approach to determine whether or not C. sp. NAI and C. anthracinus, as 
well as C. bifurcatus (gr. 1) and C. bifurcatus (gr. 2), are distinct species. However, even when using specific base 
differences, we are still faced with the same challenge inherent to other species-delimitating methods; that is, how 
much of a difference is needed for species to be considered different. 

DNA barcoding is more precise than PCR-RFLP because it allows the exact determination of base pair 
differences between individuals. Nevertheless, when a large number of individuals need to be identified, PCR- 
RFLP has been advocated as a cost-effective technique to assess molecular variation (Pfrender et al. 2010). 
However, as the cost of sequencing continues to fall, sequencing is becoming the most effective and economical 
approach, even for determining large numbers of individuals. The disadvantage of the PCR-RFLP approach is that 
if the right enzymes are not chosen for analysis, sequence nucleotide differences can go unnoticed. In our study, the 
first chosen restriction enzymes (Ssp I. Hinfl, Rsal and Taql) were not able to discern differences in the coxl 
nucleotide sequences of C. bifurcatus, C. dilutus and C. nr. atroviridis (sp. 2i) as well as those of C. staegeri, C. 
‘tigris ’ and C. frommeri. Additional restriction enzymes were necessary to separate these species. Likewise, PCR- 
RFLP did not discriminate between C. entis and C. plumosus and it would not likely be able to separate C. dilutus 
from C. pallidivittatus because both of these species pairs share identical coxl sequences. 

Overall, for non-cytological experts, we recommend the use of combined genetic and morphological 
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techniques to identify Chironomus larvae to species since this combination was much more effective than either of 
these techniques alone. 


Overall conclusions and recommendations for identifying Chironomus species 

Overall, using morphology, cytology and genetics we conclude that our 404 Chironomus larvae represent 17 
species, 14 of which have been identified as C. ( Chaetolabis ) nr. atroviridis (sp. 2i), C. ( Chaetolabis) ochreatus, C. 
{Chironomus) anthracinus, C. {Chironomus) bifurcatus, C. {Chironomus) cucini, C. {Chironomus) decorus -group 
sp. 2, C. {Chironomus) dilutus, C. {Chironomus) ends, C. {Chironomus) frommeri, C. {Chironomus) harpi, C. 
{Chironomus) maturus, C. {Chironomus) plumosus, C. {Chironomus) staegeri and C. {Chironomus) ‘tigris ’ while 
the identification of three others remains uncertain (C. sp. NAI-III). The species status of C. sp. NAI requires 
further investigation and additional studies are necessary to determine whether are not C. bifurcatus is a single 
species or a complex of at least two closely related species. Of the 14 identified Chironomus species, two belong to 
the subgenus Chaetolabis whereas 12 belong to the subgenus Chironomus. 

We collected and identified 11 (C. anthracinus, C. bifurcatus, C. cucini, C. decorus -group sp. 2, C. dilutus, C. 
entis, C. maturus, C. nr. atroviridis (sp. 2i), C. plumosus, C. staegeri and C. ‘tigris ') of the 20 Chironomus species 
currently known from the Canadian provinces east of the Rocky Mountains (that is, from Canada excluding British 
Columbia and the three northern territories; Martin 2013). Since all but three of the 31 water bodies that we 
sampled are located in the same ecozone (the Boreal Shield), some of the nine species that we did not find could be 
restricted to other ecozones, such as the prairies, where water chemistry and other factors are likely to differ from 
those in Boreal Shield lakes. In fact, it is surprising that we were able to collected so many Chironomus species 
from a single ecozone in which lake waters are generally nutrient poor, circum-neutral and soft (low concentrations 
of calcium and magnesium), which is likely to limit the range of habitats available for Chironomus species. This 
large proportion of known species is likely explained by the fact that our study lakes in this ecozone encompass 
wide ranges in these variables because some of them have been altered by discharges from mining, milling and 
smelting operations, or sewage treatment plants, or by the addition of lime to counter lake acidification (Lakes 
Arnoux, Osisko, Pelletier, Rouyn and Kelly). 

The range of chemical conditions under which some of the Chironomus species were collected was quite wide. 
For example, C. anthracinus and C. bifurcatus were found in waters that were highly acidic to circum-neutral and 
C. entis, C. bifurcatus, C. plumosus and C. staegeri were found in water bodies that were oligotrophic to eutrophic. 
In contrast, C. harpi was restricted to a highly acidic lake and C. dilutus was collected only in a lake that had been 
organically enriched by sewage. 

We found three Chironomus species (C. frommeri, C. harpi and C. ochreatus) that were previously known 
from the Nearctic (Martin 2013), but had not been reported from eastern Canada. The identification of another 
three species remains unclear (C. sp. NAI-III). We note that other species are likely to exist in eastern Canada, 
since 19 cytologically-defined but as yet unidentified or unassigned Chironomus species have been reported from 
this region (Martin 2013). Applying the combination of morphological and genetic techniques used in our study 
would likely resolve many of these taxonomic gaps in the Canadian and Nearctic fauna. 
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PROULX ETAL. 


FA07-BA10 bathophilus MC 2007 x x C. bifurcatus ( gr. 2) 

SJ07-BA1 bathophilus SJ 2007 KF278312 x C. bifurcatus (gr. 2) 

SJ07-BA2 bathophilus SJ 2007 KF278313 KF278442 x C. bifurcatus (gr. 2) 

SJ07-BA3 bathophilus SJ 2007 KF278314 KF278443 x C. bifurcatus (gr. 2) 
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FA07-SA6 salinarius MC 2007 x x C. sp. NAIII 

FA07-SA7 salinarius MC 2007 x x C. sp. NAIII 

FA07-SA8 salinarius MC 2007 x x C. sp. NAIII 

FA07-SA9 salinarius MC 2007 x x C. sp. NAIII 

. continued on the next page 








a. 

c/2 

DO 

o 

U 


DO 

£ 

o 

-C 

Cl 

i—' 

o 


=tt 

q £ 
: .2 
; c/5 

C/5 

a> 


2 ^ 
— 1 c/3 

cr 3 

D 

C/3 -g 

CD 

S <9 

b« 5 

O 


c a> 

<L> C/3 

cr 03 

<U jj> 

M *3 

’"-i CD 

^ g 
o 


a 

a 

a 

a 

a 

a 

a 

K 

a 

a 

a 

a 

a 

H 

a 

a 







< 

z 

< 

z 

< 

z 

< 

z 

< 

z 

< 

z 

<C 

z 

< 

z 

< 

z 

< 

z 

< 

z 

< 

z 

< 

z 

§ 

Z 

< 

z 

< 

z 

k 

Do 

Do 

Ik 

Do 

<5J 

k 

bo 

so 

Do 

k 

03 

Do 

03 

d. 

C/3 

d- 

C/3 

d. 

C/3 

ci 

c/3 

d- 

C/3 

d. 

C/3 

d- 

C/3 

d. 

C/3 

d. 

c/3 

d. 

CO 

d- 

C/3 

d. 

CO 

d- 

co 

d- 

co 

d. 

C/3 

d- 

CO 

>3 

,g 

to 

,g 

Co 

,g 

Co 

to 

,g 

<o 

u 

u 

u 

u 


Cj 

u 

u 

cj 

u 

O' 

0 

o 

0 

0 

o 

0 

o 

O 

0 

O 

0 


X X X X 


XXX 


X 


X 


X X 


X X X X X X X 


04 

■3’ 

— 

04 

04 

04 

3" 

3’ 

3" 

00 

oo 

OO 

O- 

o» 

o- 

OM 

04 

04 

Ph 

Ph 

a 

2 

a 

a 


O 

04 

CO 

CO 

00 

OS 

GO 

04 

o 

O' 

SO 

O' 

O' 

O' 

OO 

O' 

O- 

O 

OO 

OO 

GO 

GO 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

00 

OO 

OO 

OO 

OO 

OO 

OO 

OO 

OO 

OO 

OO 

O' 

O' 

o 

O' 

O' 

O 

O' 

O' 

o 

O' 

O' 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

04 

Ph 

Ph 

Ph 

Ph 

Ph 

Ph 

Ph 

Ph 

Ph 

Ph 

Ph 

a 

2 

a 

a 

a 

a 

a 

a 

a 

a 

2 


X X X X X 


X X X X 


X X X X 



&3 e§ e§ o§ e§ 


a: a; a; 


2 2 2 2 2 

a: a: ai cC cc 


2 £* 

<X U 


< < < < < 
Q O Q Q O 


C- 

K 


§ 

> 


3 

.3 

.3 

3 

.3 

.3 

.3 

.3 

.3 

.3 

.3 

.3 

.3 

.3 

x- 

3 

CD 

3 

c3 

CD 

CD 

CD 

n 

3 

i— 

3 

CD 

CD 

CD 

CD 

CD 

G 

G 

g 

G 

_C 

G 

G 

G 

_ 3 

c 

G 

G 

c 

_C 


CD 

.3 


*—• co 3- io *0 r" 

<<<<<< 
CO CO CO CO GO CO 



04 

CO 

3- 


SO 

o 

00 

os 

—t 

< 

< 

< 

< 

< 

< 

<C 

< 

< 

< 

C/2 

GO 


in 

C/2 



CO 

C/3 

C/2 

O' 

O 

O' 

O' 

O 

O' 

o 

O' 

O' 

O' 


^3 

r B- 


a 

CP 


s s 

3 J3 

”H. "B- 


1 1 

o- o- 


04 co 

a a 


•D- 

H-l 

CP 


CP 


3 3 ^ 


os 


< < 
« K 


< < 

x os 


< 

os 


< <C 

K K 


< 

a 


< 

(X 


g ^ 3 


J j 
< < 


IDENTIFYING CHIRONOMUS SPECIES 


Zootaxa 3741 (4) © 2013 Magnolia Press ■ 451 


DU06-PL2 plumosus DP 2006 x x C. staegeri 

DU07-PL3 plumosus DP 2007 x x C. staegeri 

DU07-PL4 plumosus DP 2007 x KF278258 x C. staegeri 

DU07-PL5 plumosus DP 2007 x x C. staegeri 
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Til I-PL2 plumosus TI 2011 KF278251 KF278399 x C. ‘tigris’ 







TABLE S2. List of voucher Chironomus species sequenced. 
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TABLE S3. Mean and range of intraspecific sequence divergences of the coxl and gbl/i genes for collected and reference Chironomns species. 
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TABLE S5. Bases that differ between coxl sequences of C. bifurcatus groups 1 and 2. Refer to Fig. 8 for sequence label. 


Base position 


Species 

Sequences 

1 

2 

3 

3 

3 

4 

4 

4 

5 

6 



9 

7 

3 

4 

8 

3 

7 

9 

1 

2 



9 

1 

7 

3 

2 

3 

8 

9 

7 

2 

C. bifurcatus (gr. 1) 

bathophilus (n=l / DA / KF278315) 

C 

C 

T 

T 

A 

G 

T 

C 

C 

A 

C. bifurcatus (gr. 1) 

bathophilus (n=l 1 / AD, AR, DP, OP / KF278316- 
KF278326) 

C 

C 

T 

T 

A 

G 

T 

C 

C 

A 

C. bifurcatus (gr. 1) 

C. bifurcatus (KF278361) 

c 

c 

T 

T 

A 

G 

T 

A 

c 

A 

C. bifurcatus (gr. 1) 

bathophilus (n=l / SJ / KF278352) 

c 

c 

T 

T 

A 

G 

T 

C 

c 

A 

C. bifurcatus (gr. 2) 

C. bifurcatus (KF278353) 

T 

A 

A 

C 

T 

A 

A 

T 

T 

C 

C. bifurcatus (gr. 2) 

C. bifurcatus (KF278345) 

T 

A 

A 

C 

T 

A 

A 

T 

T 

C 

C. bifurcatus (gr. 2) 

bathophilus (n=l / KI / KF278307) 

T 

G 

A 

c 

T 

A 

A 

T 

T 

c 

C. bifurcatus (gr. 2) 

bathophilus (n=3 / SJ / KF278312-KF278314) 

T 

G 

A 

c 

T 

A 

A 

T 

T 

c 

C. bifurcatus (gr. 2) 

bathophilus (n=2 / MC / KF278310-KF278311) 

T 

G 

A 

c 

T 

A 

A 

T 

T 

T 

C. bifurcatus (gr. 2) 

bathophilus (n=2 / MC, TI / KF278308-KF278309) 

T 

G 

A 

c 
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A 

A 

T 

T 
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